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(Dated: February 2, 2009) 

Abstract 

Using ^(3770) DD events collected with the CLEO-c detector at the Cornell e~^e~ storage 
ring, tagged by fully reconstructing one D meson in a hadronic decay mode, we measure absolute 
branching fractions and differential decay rates for — > 7r~e+fe) —>■ Tr'^e+fg, 
and Z)+ K^e^Ve- The measured decay rates are used to study semileptonic form factors gov- 
erning these transitions and to test unquenched Lattice QCD (LQCD) calculations. We average 
our results with previously published CLEO-c measurements of the same quantities using a neu- 
trino reconstruction technique. Combining LQCD calculations of form factor absolute normaliza- 
tions /+(0) and measurements of /^(0)|V;d| and (O)!^^!, we find iKdl = 0.222(8)(3)(23) and 
I Fes I = 1.018(10) (8) (106), where the uncertainties are statistical, experimental systematic, and 
from LQCD, respectively. 
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I. INTRODUCTION 



The quark mixing parameters are fundamental constants of the Standard Model of parti- 
cle physics. They determine the nine weak-current quark coupling elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [1]. In the Standard Model the CKM matrix is unitary. 
Measuring the quark couplings tests the unitarity of the matrix. 

The extraction of the quark couplings is difficult because quarks are bound inside hadrons 
by the strong interaction. Semileptonic decays are the preferred way to determine the CKM 
matrix elements as the strong interaction binding effects are confined to the hadronic current. 
They are parameterized by form factors that are calculable, for example, by lattice quan- 
tum chromodynamics (LQCD) and QCD sum rules. Nevertheless, form factor uncertainties 
dominate the precision with which the CKM matrix elements can be determined [2]. 

Studies of the semileptonic decays of D mesons play an important role in understanding 
the CKM matrix. First these decays allow the robust determination of the couplings iV^^j 
and \Vcd\ by combining measured branching fractions with form factor calculations. Second 
\Vcs\ and \Vcd\ are tightly constrained when the CKM matrix is assumed to be unitary. 
Therefore measurements of charm semileptonic decay rates, when combined with the values 
of \Vcs\ and \ Vcd\ constrained by the unitarity of the CKM matrix, rigorously test theoretical 
predictions of D meson semileptonic form factors. 

Recently using ■?/'(3770) DD events and a neutrino reconstruction technique combined 
with an independent measurement of the number of D mesons, CLEO reported the most 
precise determinations of the absolute branching fractions and differential decay rates dV/dq^ 
for the decays n'e+v^, D+ 7r°e+z/e, i^"e+z/e, and D+ -> K^e+Ve [3]. 

(Throughout this paper charge- conjugate modes are implied.) The differential decay rates 
were used to determine the absolute magnitude and shape of the semileptonic form factors 
and to determine \Vcs\ and \Vcd\- In this paper we present a complementary analysis which 
measures the same quantities with similar precision in a common data set but with a different 
technique that is independent of the number of D mesons in the data sample. The two 
analyses obtain consistent results, providing increased confidence in their correctness, and 
each represents a marked improvement in our understanding of charm semileptonic decays. 

As the two analyses use a common data set, the results are correlated. We calculate 
average values of the branching fractions, form factors and \Vcs\ and \Vcd\ measured in the 
two analyses, taking into account correlations between them. The average values represent 
the best determinations of these quantities with the CLEO-c 281 pb~^ data set. 

The paper is organized as follows. We review the semileptonic decay formalism in Sec. II. 
The data sample and CLEO-c detector are described in Sec. III. The analysis technique to 
identify semileptonic decays is introduced in Sec. IV. In Sees. V and VI we describe the use 
of this technique to measure the absolute branching fractions, differential decay rates and 
form factor parameters for {D^) decays to 7r~e^z/e {^K^e^Uf.) and K^e^v^ {K^e^v^. The 
extraction of CKM parameters is described in Sec. VII. In Sec. VIII we average the results 
presented here with the results obtained in [3] . Finally, in Sec. IX a summary is provided. 

II. SEMILEPTONIC DECAY FORMALISM 

The matrix element for a semileptonic decay Mj(gjg') — >■ Mf[qfq')i^i'£ where Mj and Mf 
are the initial and final state mesons, and qf are the initial and final state quarks, and q' 
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is a spectator anti-quark, can be written as 



M{M, ^ Mf£+ue) = -&V*L^^H„ (1) 



V2 



where Gi? is the Fermi constant, Vq^q^ is the appropriate CKM matrix element, and L'^ and 
Hfj^ are the leptonic and hadronic currents, respectively. The leptonic current is known and 
can be written in terms of the lepton and neutrino Dirac spinors, ue and v^,, 

L^^ = ua^{l-^,)v,. (2) 

The underlying simplicity of the weak transition qi qjW^ is obscured by the strong 

interaction as the initial and final state quarks are bound within hadrons. The hadronic 
current can be written as 

H,^{M^\qf^^{l-^^)q^M,). (3) 

The hadronic current describes the non-perturbative strong interaction physics of hadron 
formation. Usually, one exploits the fact that the hadronic current transforms as a four vector 

under Lorentz transformations by parameterizing it with a set of invariant form factors. 
This is achieved by constructing all possible quantities with transformation properties of 
four vectors from the momenta of particles involved in the decay, their spin - polarization 
vectors and invariant tensors, and expanding the hadronic current in terms of these with an 
invariant form factor multiplying each of them. The form factors can only be functions of 
Lorentz scalars. In Mj — * Mfi'^u^, there is one such invariant, which is usually chosen to be 
q"^, the square of the invariant mass of the virtual W. 

In pseudoscalar-to-pseudoscalar semileptonic decays {Pi{qiq) — > Pf{qfq)£~^i'i), the 
hadronic current has a simple structure: 



{Pf{pf)\V^\PM) 

mt — m 



+Mq')-^^^{p^-Pfr, (4) 



where pi {rrii) and pf {mf) are the four-momenta (masses) of the initial Pi and final Pf 
mesons, and /+(?^) and fo{Q'^) are the form factors governing the transition. Kinematic 
constraints require /+(0) = /o(0). In the limit of negligible lepton mass, which is applicable 
for £ — e, only one form factor remains, 

{PfiPfW'^lPM) = fWKp^ + PfT. (5) 

The form factor f+iq"^) measures the probability to form the final state hadron; it is largest 
when the daughter meson is stationary in the parent meson rest frame = g^ax) 
smallest when the daughter meson is moving with maximum velocity in the parent meson 
rest frame q'^ — 0. 

The differential decay rate is given by 

dr _ Gp\Vq^qj^\'^Pp 



dq^ 247r3 



i/+(or, (6) 
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where pp^ is the magnitude of the three-momentum of the Pf meson in the rest frame of 
Pi. The shape of the distribution is dominated by the dependence on p|>^, which arises 
because the decay proceeds via a P-wave. This dependence significantly enhances the rate 
at low g^. We perform fits to the differential decay rate to measure the four semileptonic 
modes D Ke^v and D ne^u. In this paper we denote the form factor governing 
D Ke^u and D vre+i/ by f^{q^) and f^{q'^), respectively. 



A. Parametrization of the Form Factor Dependence 



The dependence of the form factors on is unknown, as it is determined by non- 
perturbative QCD. One may express the form factors in terms of a dispersion relation, 
an approach that has been well established in the literature (see, for example, Ref. [4] and 
references therein): 

u u) = ^ ' , + i r ""'^:"».it (7) 

■'^y^J 1 - A . q^ -K Jimn+mpYt-q^ -is ^' 

where Mpoic is the mass of the lowest lying (qiqf) meson with the appropriate quantum 
numbers: for D Ke^u^ it is D*'^{1~) and for D Tie^Ve it is D*+(l^), the parameter A 
gives the contribution from the vector pole at = 0, mD is the mass of the D meson, and 
mp is the mass of the final state pseudoscalar meson. It is common to write the dispersive 
representation in terms of an explicit pole and a sum of effective poles, 

u,,^).m^-±^.j:^-, (8, 

^ ~ 1 %— k=l 1 2.^ 

where pk and 7^ are expansion parameters that are not predicted. 

A series expansion around (f — to [5-8], where to is defined below, is commensurate with 
the dispersion relations. As expansions in suffer from convergence problems due to the 
presence of nearby poles, the expansion is formulated as an analytic continuation into the 
t = q^ complex plane. There is a branch cut on the real axis for t > M^,^ corresponding to 
D{K, tt) production, that is mapped onto the unit circle by the variable z defined as 



z{q^,h) 



^± = {niD^mpf, (9) 



where to is the arbitrary g value that maps to ^ = 0. We choose to = — y 1 — t-/t+) 
because this choice minimizes the maximum value of z in the decay (|-Smax| = 0.051 for D 
Ke'^Ve and |2;niax| = 0.17 for D — >• Tie'^Ve)- 
The form factor is given by 

= . . fl + E auitMe. t'f] , (10) 



P{q')HQMo) 



k=l 
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where P(g^) = 1 for D — tt and Piq^) — z{(f ,m'i^»+) for {D — > K), and is arbitrary. 
Physically P accounts for the presence of the pole, and is chosen to enable a simple 
expression for the series in terms of the a^. We follow Ref. [8]: 



n<i ,to) = c' 




(11) 



This choice leads to the constraint 



E«^<1 (12) 

k=l 

for any choice of ric- To leading order the coefficient c is given by 



7iml/3, (13) 



where is the charm quark mass, which we take to be 1.2 GeV/c . An advantage of 
the z expansion is that it is model independent and satisfies analyticity and unitarity. In 
addition, measuring the in D — > iri'^i'i constrains the class of form factors needed to fit 
B — > iri'^Uj^ and hence may improve the determination of \Vub\. Finally, in Heavy Quark 
Effective Theory (HQET) [9] there exist relations between the Oj in D and B semileptonic 
decays. 

The expansion parameters are not predicted. As z is small, the series is expected to 
converge quickly. Recently BABAR [10], using a data sample of 75,000 D° — > K~e'^Ug 
events, found the differential rate to be well described with only a linear term. In this work 
we will fit the data to both linear and quadratic terms and use the series expansion for our 
main results. There are alternatives to the z expansion [11]. 

In order to compare to lattice QCD calculations and previous measurements, we will also 
compare the data to other parametrizations of the form factor dependence. A variety of 
models have been traditionally used to parameterize the dependence. The most common, 
based on vector meson dominance [4], uses only the first term in the dispersion relation. In 
this "simple pole model" the q^ dependence is given by 

U{g') = . (14) 

pole 

Previous measurements of the spectrum in K^i^Ui, the best measured charm 

semileptonic decay, find a value of the pole mass many standard deviations from Md* [3, 
10, 12-14]. At low to medium values of q"^ the q'^ spectrum is distorted compared to a simple 
pole, suggesting contributions from a spectrum of poles above the pole with the lowest mass. 

The modified pole or Becircvic-Kaidalov (BK) parametrization [15] attempts to address 
the shortcoming of the simple pole model by keeping the first term in the dispersion relation 
sum. The form factor is given by 

) = -J^ (16) 

pole pole 
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where Mpoie is the spectroscopic pole mass and a, a free parameter, is an additional "effec- 
tive" pole which represents the total contribution of all additional poles. 

In current data the evolution of form factors are indistinguishable from straight lines. 
Therefore it is convenient to define the physical shape observables in terms of form factor 
slopes at i = [16, 17] 



1 

5 = 




t=o 



dt 



(16) 



t=o> 



The quantities /3 and 5 depend on the masses of the mesons involved, and as they arc physical 
quantities they are independent of the renormalization scale or scheme. 

The BK parametrization requires several assumptions to reduce the multiple parameters 
initially present (Eq. (8)) to one. Specifically, it is assumed that f3, which measures scaling 
violations, is near unity, and 6, which measures spectator quark interactions, is near zero. 
This sets the physical observable 



1 + 1/P-S 



m 



D 



m, 



/+(0) dq^ 



2 (at q^ = 0), 



(17) 



as noted in Rcf. [17], corresponding to a ~ 1.75 for D Ki'^i'i and 1.34 for D tiC^v^. 
Previous experimental measurements of the q^ spectrum in D — > Kj'nl'^Vji do not agree with 
this value of a [3, 10, 12-14]. 

Although the simple pole model and modified pole model are unable to describe the 
spectrum of the data when the pole mass is fixed to the relevant spectroscopic pole, or 
OL ~ 1.75 for D Kl^Vf and 1.34 for D — > nl^vp^ they do describe the data well for values 
of the shape parameters many standard deviations from the expected values. 



B. Form Factor Calculations 

A variety of model dependent calculations of form factors exist. In these models the form 
factors are evaluated at a fixed value of q^, e.g., (f — or (f — q^^ = {niD — tupY, and 
are extrapolated over the full range of q^ using a parametrization, such as those discussed 
above. 

Quark model calculations estimate meson wave functions and use them to compute the 
matrix elements that appear in the hadronic current. There arc a large variety of theoretical 
calculations [18]. Among them the ISGW model [19] has been widely used to simulate heavy 
hadron semileptonic decays. This model is expected to be valid in the vicinity of q^ = q^^, 
the region of maximum overlap between the initial and final meson wave functions. 

In the ISGW model the form factors are assumed to have the form 

f{<f) = fiQlJe-'^^'---'^^- (18) 

The ISGW2 model [20], an update of the ISGW model, incorporates constraints from heavy 
quark symmetry. It uses a dipole term for the form factor q^ dependence expressed in terms 
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of the radius of a meson (r) rather than the mass of the appropriate (qiqf) meson: 



Mq') = MqU (l + ^ (da. - q')) ' ■ (19) 

The ISGW2 model predicts (g^ax) = 1-23 and = 1.12 GeV"^ [20]. Previous measure- 
ments (e.g., Refs. [10, 12-14]) do not agree with these values. 

QCD sum rules [21, 22], are expected to be valid at low q^. For — > K~i~^i'e, and 
using a value of 150 MeV for the strange quark mass, one obtains [22] /:^(0) = 0.78(11) and 
aK = 0.07t{|:Jf using the modified pole ansatz. For D ni+Ui [22] reports /^(O) = 0.65(11) 
and = O.Oltg:^^. 

The above models are based on theoretical assumptions and, in consequence, introduce 
a difficult to quantify theoretical uncertainty that is significantly larger than the presently 
achievable experimental statistical and systematic uncertainties combined. Therefore this 
limits the precision with which \ Vcs\ and \ Vcd\ can be determined from exclusive semileptonic 

charm meson decays. 

Lattice QCD computes f+{q^) from first principles. Current results must be extrapolated 
to physical values of light quark masses and corrected for finite lattice size and discretization 
effects. There have been several evaluations of f+{q'^) for different values of the momentum 
transfer in the quenched approximation [23, 24]. These results, which do not include QCD 
vacuum polarization, have been combined [23], to give f^{0) = 0.73(7). LQCD calculations 
which incorporate QCD vacuum polarization (unquenched calculations) have produced re- 
sults that agree with experiment to within a few percent for a number of quantities [25]. 
The first unquenched LQCD calculation [26] of form factors in D ^ Ke~^u,. and D Tie^Ue 
reports f^{Q) = 0.73(3)(7), qk = 0.50(4), /^(O) = 0.64(3)(6), and = 0.44(4) using 
the modified pole ansatz to parameterize the q^ dependence of the form factor. Here the 
systematic uncertainty is dominated by the effect of discretization. While the form factors 
are currently calculated to a modest precision of ten percent, the uncertainties are system- 
atically improvable to a precision that matches, or exceeds, the experimental measurements 
presented here and in [3]. Accordingly, we use [26] to extract values for |Vcs| and \Vcd\ in 
this work. 

III. DATA SAMPLE AND THE CLEO-c DETECTOR 

The data sample used in this analysis consists of 281 pb~^ of e~^e~ annihilation data 
taken at the ■?/'(3770), which is about 40 MeV above the DD pair production threshold. 
(Throughout this paper D is used to denote and D^.) The data include approximately 
1.0 X 10^ L>°L>° events and 0.8 x 10^ D+D- events. 

CLEO-c is a general-purpose solenoidal detector. The charged particle tracking system 
covers a solid angle of 93% of An and consists of a small-radiiis six-layer low mass stereo wire 
drift chamber concentric with and surrounded by a 47-layer cylindrical drift chamber. The 
chambers operate in a 1.0 T magnetic field and achieve a momentum resolution of ~0.6% at 
p —1 GeV/c. The main drift chamber provides specific-ionization [dE/dx] measurements 
that discriminate between charged pions and kaons. Additional hadron identification is 
provided by a Ring-Imaging Chcrcnkov (RICH) detector covering approximately 80% of 
47r. Identification of positrons and detection of neutral pions rely on an electromagnetic 
calorimeter consisting of 7800 cesium iodide crystals and covering 95% of 47r. The calorimeter 
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achieves a photon energy resohition of 2.2% at —1 GeV and 5% at 100 MeV. The 
CLEO-c detector is described in detail elsewhere [27]. 

The response of the CLEO-c detector was studied using a GEANT-based [28] Monte 
Carlo (MC) simulation. To develop selection criteria and test the analysis technique several 
MC simulations are used. '?/'(3770) — > DD events are generated using EvtCen [29] and 
each D meson is allowed to decay in accordance with the best experimental and theoretical 
information. We refer to this as "generic MC" . The MC sample generated corresponds to an 
integrated luminosity of about 11 fb~^ which is a factor 40 larger than the data. Semileptonic 
signal decays are generated with the modified pole model form factors [15] with parameters 
from the most recent unquenched LQCD calculations [26]. 

Due to the tagging technique employed in the analysis, backgrounds from the non-DD 
processes e+e^ — > qq, where q is a. u, d, or s quark, e"'"e~ — >• r+r", and e"'"e~ — >• ip{2S)'y, 
are nearly absent. These non-DD processes are also modeled using MC simulation and are 
scaled absolutely according to their measured cross sections at the ■?/'(3770). 

A second type of MC sample, which we refer to as "signal MC" , consists of several samples 
of '0(3770) — > DD events in which the D is allowed to decay to all possible final states, and 
the D decays to a specific semileptonic final state. 



IV. EVENT RECONSTRUCTION 

The reconstruction technique used in this analysis was first applied by the Mark III 
collaboration [30] at SPEAR. This technique was used to measure D semileptonic branching 
fractions with a smaller data sample at CLEO-c [31]. That data sample was too small to 
study charm semileptonic form factors, which are the focus of studies reported in this paper. 

The presence of two D mesons in a ^(3770) event allows a tag sample to be defined in 
which a D is reconstructed in a hadronic decay mode. A sub-sample is then defined in 
which a positron and a set of hadrons, as a signature of a semileptonic decay, are required in 
addition to the tag. Tagging a D meson in a -0(3770) decay provides a D with known four- 
momentum, allowing a semileptonic decay to be reconstructed with no kinematic ambiguity, 
even though the neutrino is undetected. 

The tag yield can be expressed as iVtag = ^Nd D^tsig^tag, where Ndd is the produced 
number of DD pairs, i3tag is the branching fraction of hadronic modes used in the tag 
sample, and ctag is the tag efficiency. The yield of tags with a semileptonic decay can be 
expressed as iVtag,sL = ^NDoBtagBsL^tag.sL where Bsl is the semileptonic decay branching 
fraction, including subsidiary branching fractions, and etag,sL is the efficiency of finding the 
tag and the semileptonic decay in the same event. From the expressions for N^g^g and A^tag,SL 
we obtain 

^tag,SL £tag -^tag,SL/ £ /^^x 

^SL = -77 ; = —77 , (20) 

where e = etag,SL/etag is the effective signal efficiency. The branching fraction determined 
by tagging is an absolute measurement. It is independent of the integrated luminosity and 
number of D mesons in the data sample. Due to the large solid angle acceptance and high 
segmentation of the CLEO-c detector and the low multiplicity of the events etag,SL ~ ctagCsL, 
where esL is the semileptonic decay efficiency. Hence the ratio etag,sL/etag is insensitive to 
most systematic effects associated with the tag mode and the absolute branching fraction de- 
termined with this procedure is nearly independent of the tag mode. Below, we first describe 
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the procedure used for the reconstruction of tags followed by that for the reconstruction of 
semileptonic decays [31]. 

A. Tag Selection 

Hadronic tracks must have momenta above 50 MeV/c and | cos6'| < 0.93, where 9 is the 
angle between the track direction and the beam axis. Identification of hadrons is based 
on measurements of specific ionization in the main drift chamber and information from the 
RICH. Pion and kaon candidates are required to have dE/dx measurements within three 
standard deviations (3cr) of the expected value. For tracks with momenta greater than 
700 MeV/c, RICH information, if available, is combined with dE/dx. The efficiencies (95% 
or higher) and misidentification rates (a few percent) are determined with charged pion and 
kaon samples from hadronic D decays. 

We select tt" candidates from pairs of photons, each having an energy of at least 30 MeV, 
and a shower shape consistent with that expected for a photon. A kinematic fit is performed 
constraining the invariant mass of the photon pair to the known 7r° mass. The candidate is 
accepted if the unconstrained invariant mass is within 3(J, where a (typically 6 MeV/c^) is 
determined for that candidate from the kinematic fit, and the kinematic parameters for the 
7r° determined with the fit are used in further reconstruction. 

Candidate events are selected by reconstructing a or D~ tag in the following hadronic 
final states: K+n', K+n-ix^, K+n-n^n^, K+ti-'k-'k+, ^sTt", ir^7r-7r+, K^tt-tt+ttO, and 
K~ for neutral tags, and /T^tt", K^tt~ti~^ Kgir'n'^, K~^iT~7r~7r^ , K^tt~tt~ti^^ and 
K~ K^n" for charged tags. These modes constitute about 46% and 28% of all and 
D~ decays, respectively. Tagged events are selected using two variables: /S.E = Ed — -Ebeam, 
the difference between the energy of the tag candidate {Ed) and the beam energy (-Bbeam), 

and the beam-constrained mass Mbc = \/ E^cmn/ '^'^ " \Pd\'^/c^, where is the measured 
momentum of the tag candidate. Note that the use of -Ebcam instead of Ed improves the 
resolution of Mbc by one order of magnitude, to about 2 MeV/c^, which is dominated by 
the beam energy spread. If multiple candidates are present in the same tag mode, the one 
candidate per tag flavor with the smallest AE is chosen. 

The number of tags reconstructed in each mode is obtained by imposing a mode dependent 
requirement on AE, counting the number of events in the signal region of Mbc, defined as 
-6.5 MeV/c^ < (Mbc-mip) < 9.5 MeV/c^, where uid [32] is the known D meson mass, and 
subtracting the background contribution from it. Fits to the Mbc distributions, shown in 
Figs. 1 and 2, are made using the procedure described in [33]. We fit the Mbc distributions 
to a signal shape and one or more background components. The signal shape includes the 
effects of beam energy smearing, initial state radiation, the line shape of the ip{3770), and 
reconstruction resolution. The background is described by an ARGUS function [34] , which 
models combinatorial contributions. The background contribution in the signal region is 
estimated by integrating this function. The yields of the eight neutral tag modes and the six 
charged tag modes, and their reconstruction efficiencies as determined with the generic MC 
simulation, are given in Tables I and II. There are approximately 3.1 x 10^ and 1.6 x 10^ 
neutral and charged tags, respectively. 
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FIG. 1: Fits (solid line) to the Mbc distributions in data for eight tag modes. The backgrounds 
are shown by the dashed line. 



TABLE I: Yields with statistical uncertainties and reconstruction efficiencies of D tags. 



Tag 


; Mode 




^tag (%) 




-> K+TT- 


51002(230) 


64.77(3) 






98117(347) 


33.30(1) 






23040(220) 


14.41(1) 


L>o 




77641(303) 


45.46(2) 


^" 




24533(187) 


38.33(2) 






20355(260) 


17.81(5) 






8175(99) 


31.01(5) 




K-K+ 


4614(76) 


57.35(9) 


All Neutral Tags 


307478(657) 
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FIG. 2: Fits (solid line) to the distributions in data for six D tag modes. The backgrounds 
are shown by the dashed line. 

TABLE II: Yields with statistical uncertainties and reconstruction efficiencies of D~ tags . 



Tag Mode iVtag etag (%) 



D- 


— > K'^n TT 


79896(291) 


53.81(2) 


D- 


^ K+TT-ir-iT^ 


23740(196) 


25.23(2) 


D- 




11456(113) 


45.14(5) 


D- 




25159(210) 


21.97(2) 


D- 




16431(191) 


31.58(3) 


D- 


K-K+TT- 


6794(100) 


44.72(5) 



All Charged Tags 163476(477) 



B. Selection of Semileptonic Decays 

After a tag is identified, we search for a positron and a set of hadrons recoiling against the 
tag. (Muons are not used as D semileptonic decays at the V'(3770) produce low momentum 
leptons for which the CLEO-c muon identification is not efficient.) Positron candidates are 
required to have momenta of at least 200 MeV/c and to satisfy | cos^^| < 0.90, where 9 is 
the angle between the positron direction and the beam axis. The efficiency for positron 
identification rises from about 50% at 200 MeV/c to 95% just above 300 MeV/c and is 
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roughly constant thereafter. The rate for misidentifying charged pions and kaons as positrons 
averaged over the momentum range is approximately 0.1%. The energy lost by positrons to 
bremsstrahlung photons is partially recovered by adding showers that are within 5° of the 
positron momentum and are not matched to other particles. The selection of 7r~, 7r°, K~ , 
and candidates is identical to that used for tags. 

The tag and the semileptonic candidate are then combined. Events that include tracks 
other than those of the tag and the semileptonic candidate are vetoed [35]. After all selection 
criteria are applied, multiple candidates in the same event are rare in all modes except 
— > TT^e+i^e. For — > 7r°e+z/e, in the few percent of events with multiple candidates, 
one combination is chosen per tag candidate based on the proximity of the invariant masses 
of the 7r° candidates to the expected mass. 

Semileptonic decays are identified using the variable U = i^miss ~ clpmissl, where ^^miss 
and Pmiss ^-re the missing energy and momentum of the D meson decaying semileptonically, 
calculated using the difi^crcncc of the four-momentum of the tag and that of the observed 
products of the semileptonic decay. If the decay products of the semileptonic decay have 
been correctly identified, U is expected to be zero, since only a neutrino is undetected. To 
improve the resolution in U, the crossing angle of the beams (~ 3 mrad) is allowed for 
by recalculating all track momenta and shower energies in the ( (3770) rest frame, and the 

four-momentum of the tag is approximated by (-Ebeam/c, y''(£'beam/c)^ — {cmDYPD)-i where 
is the unit direction vector of the D in the '0(3770) rest frame determined using the 
direction of the D tag in the same frame. Due to the finite resolution of the detector, the 
distribution in U is approximately Gaussian, centered at = with a ^ 12 MeV, for all 
modes except D'^ Tr'^e'^Ue, for which a is approximately two times larger. 

Using this procedure we obtain the U distributions shown in Fig. 3. For each mode a 
clear signal is evident centered on U = 0, while backgrounds are very small near U = 0. 
In —>■ 7r~e+i/e the peak at positive U is from two sources: D° — > K~e'^i'e when a K~ 
is misidentified as a 7r~ (peak at 130 MeV) and from D — > K'n'^n^ where the K~ is 
mistaken for an electron and the 7r° is unobserved (peak at 180 MeV). This background is 
present because each event is not required to have both a D° and a D'^. Specifically, on the 
semileptonic side of the event both vr^e+z/g and 7i^e~Ue are accepted, on the 

tag side for example both — > i^+7r~ and D° — > K~7i^ are accepted. The kaon produced 
in the decay of the tag is not required to have the same charge as the lepton produced in the 
semileptonic decay. If this requirement were made the D K~tt^h^ background would be 
removed, but decay sequences where the tag undergoes a doubly Cabibbo suppressed decay 
such as D'^ — > K~^7T~ , and D'^ — > 7T~^e~Ue would be removed as well. 

The yield for each semileptonic mode is determined from a fit to the corresponding U 
distribution, as shown in Fig. 3 with all tag modes combined. The yields are reported in 
Table III. In each case the signal function consists of a Gaussian to describe the core of the 
U distribution and two power law tails to account for initial and final state radiation (ISR 
and FSR): 



ai(^-Q;i+i)-"i t>ai 

/ = S e-^ -a2 <t <ai , (21) 

. a2(^ -a2- t)-^' t < -a2 

where i = (C/— C/mean)/c"c/, 0,1 = (ni/ai)"^ e""i/^, and 02 = (n2/a2)"^^ e'"'^/^. The parameters 
describing the tails of the signal function (cci, 0:2, rii, and 712) are always fixed in fits to the 
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FIG. 3: Fits (solid line) to the U distributions in data (points) for D*^ — 7r~e+i/e, -D° — K~e''^Ue, 
— *■ 7r°e"''fe, and D'^ — > K^e^Ve- The background contributions are represented by dotted or 
dashed lines. In 7r~e"'"t'e the background peaks at positive U are described in the text. 

data to the values found in signal MC simulation. The au is fixed to the value predicted by 
the MC simulation in the fit for D° — > 7r°e"'"fe, which has the smallest signal yield and the 
largest background level among the four semileptonic modes, and allowed to float in the flts 
for the other modes. 

The background functions are determined from the generic MC simulation. The back- 
grounds are small and arise mostly from misrcconstructed semileptonic decays with correctly 
reconstructed tags. The background shape parameters are flxed, while the background nor- 
mahzations are allowed to float in all fits to the data. 

V. ABSOLUTE BRANCHING FRACTION MEASUREMENTS 
A. Determination of the Branching Fractions 

The absolute semileptonic branching fractions are obtained from our tagged semileptonic 
yields A^tag,SL, tag yields A^tag, and the efficiencies e, using Eq. (20). The simulation of 
each semileptonic mode employs the simple pole model with Mpoic = 2.0 GeV/c^. The 
efficiency depends weakly on Mpoic; accordingly the efficiencies are re- weighted to the value 
of Mpoie measured in the data. These efficiencies are then weighted by the tag yields shown in 
Tables I and II to obtain the overall efficiency. The absolute semileptonic branching fractions 
are obtained using these weighted efficiencies. Table III presents our absohite semileptonic 
branching fraction measurements with statistical and systematic uncertainties. A description 
of how the systematic uncertainties are obtained is provided in the next subsection. 

The procedure for measuring semileptonic branching fractions is tested using the generic 
MC sample. In the test, the MC sample is treated identically to the data. In addition, 
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the procedure was separately tested for each combination of tag and semileptonic mode. 
We find that the input and output branching fractions are consistent within statistical 
uncertainties in all cases. The largest deviation is observed for — > K'^e^Ue with all tag 
modes combined, where the discrepancy is less than one third of the statistical uncertainty 
on the measurement. 

To check the consistency of the measurement of the semileptonic branching fractions, we 
have also measured semileptonic branching fractions for each tag mode separately for the 
two Cabibbo allowed final states where there are adequate statistics in each tag mode. We 
present the results in Tables IV and V. 

We note that the effective semileptonic efficiency is larger for tag modes with higher 
multiplicity. This happens primarily because tag reconstruction efficiencies in events with 
the second D meson decaying hadronically are slightly smaller compared to signal events 
with the second D meson decaying to a low multiplicity semileptonic final state. 

We find that the branching fractions arc consistent among tag modes. The results in 
Tables IV and V also demonstrate consistency between the weighted averages of the indi- 
vidual tag mode branching fractions and the branching fractions obtained with all tag modes 
combined. 

TABLE III: Signal efficiencies, yields, and branching fractions in this work (first four columns) 
and, for comparison, the branching fractions measured using the first 56 pb'^ CLEO-c V(3770) 
data sample [31], and values from PDG-04 [32]. The first uncertainty is statistical and the second 
systematic in the fourth and fifth columns, and statistical or total in the other columns. 



Decay Mode 




-^tag.SL 


^SL (%) 


fisL (%) (56 pb-i) 


HsL (%) (PDG-04) 




72.54(11) 


699(28) 


0.314(13)(4) 


0.262(25)(8) 


0.36(6) 


D+ 7r°e+i/e 


44.72(13) 


281(19) 


0.384(27)(23) 


0.44(6)(3) 


0.31(15) 


DO ^ K~e+Ve 


61.06(7) 


6786(84) 


3.61(5)(5) 


3.44(10) (10) 


3.58(18) 


D+ K^e+Ve 


20.01(4) 


2910(55) 


8.90(17)(21) 


8.71(38)(37) 


6.7(9) 



B. Study of Systematic Uncertainties for Absolute Branching Fractions 

We have considered the following sources of systematic uncertainty in the measurements 
of branching fractions and give our estimates of their magnitudes in parentheses. The un- 
certainties associated with the efficiency for finding a track (0.3% for each pion, kaon, or 
positron, combined in quadrature with an additional 0.6% for each kaon), for reconstructing 
a (4.3%), and for reconstructing a (1.8%), are estimated using missing mass tech- 
niques described in [33]. The uncertainty in the positron identification efficiency (1.0%) 
is obtained using a comparison of the detector response to positrons from radiative pro- 
cesses in the data and MC simulation. The effect of the event complexity is incorporated 
by studying positrons both in isolation and embedded in hadronic events. Uncertainties in 
the charged pion and kaon identification efficiencies (0.1% per pion and 0.2% per kaon) are 
estimated using hadronic D meson decays. The uncertainty in the number of tags (0.4%) 
is estimated by using alternative signal functions in the fits to the distributions and 



15 



TABLE IV: D K~e~^i'e semileptonic yields in data, the semileptonic efficiency for each hadronic 
tag mode, e, and the branching fraction measurement for each hadronic tag mode. The last two 
lines show the weighted average of the individual measurements and the result from the fit with 
all tag modes combined. 



Mode 


A^tag,SL 


e (%) 


BsL (%) 




1088(34) 
2143(47) 
593(25) 


59.36(15) 
61.66(12) 
67.11(30) 


3.60(11) 
3.55(8) 
3.84(15) 




1693(42) 


59.44(13) 


3.67(9) 




516(23) 


59.47(24) 


3.54(16) 


i^g7r-7r+7r° 


474(22) 


64.46(30) 


3.61(17) 


K-K+ 


160(13) 
118(11) 


60.52(43) 
59.45(51) 


3.23(26) 
4.32(42) 


Average: 






3.61(5) 


Combined Fit: 


6786(84) 


61.06(7) 


3.61(5) 



TABLE V: K^e^Ve semileptonic yields in data, the semileptonic efficiency for each hadronic 

tag mode, e, and the branching fraction measurement for each hadronic tag mode. The semileptonic 
efficiency includes subsidiary branching fractions [36] . The last two lines show the weighted average 
of the individual measurements and the result from the fit with all tag modes combined. 



Mode 


-^tag,SL 




BsL (%) 




1437(39) 


19.88(5) 


9.04(25) 




430(21) 


20.51(9) 


8.83(44) 




201(14) 


19.91(18) 


8.81(46) 


if^TT-TrO 


443(22) 


20.17(9) 


8.73(43) 




272(17) 


19.82(11) 


8.35(53) 


K-K+TT- 


130(12) 


19.97(16) 


9.59(88) 


Average : 






8.89(17) 


Combined Fit: 


2910(55) 


20.01(4) 


8.90(17) 



by varying the end point of the background function [34]. The uncertainty associated with 
the requirement that there be no additional tracks in tagged semileptonic events (0.3%) is 
estimated by comparing fully reconstructed DD events in data and MC simulation. The 
uncertainty associated with the number of signal events is estimated by using an alternative 
signal function (a double Gaussian) in the fits and by counting events in the signal region 
(4.2% for TT^e'^i'e, 0.3% for all other modes). The uncertainty in the semileptonic 

reconstruction efficiencies due to imperfect knowledge of the semileptonic form factors (0.0% 
to 0.3% depending on mode) is estimated by varying the form factor shape parameters in 
the MC simulation within uncertainties in their measurements reported in Sec. VI D. The 
uncertainty associated with the simulation of FSR and bremsstrahlung radiation in the de- 



16 



TABLE VI: Summary of systematic uncertainties considered in the measurements of absolute 
branching fractions of the four semileptonic modes. The modes are labeled by their final state 
hadrons. 



Systematic uncertainty (%) 







TT 






Number of D tags 


0.4 


0.4 


0.4 


0.4 


Electron ID efficiency 


1 n 
i.U 


i.U 


1 n 
i.U 


i.U 


Hadron ID efficiency 


0.2 


0.1 


0.0 


0.0 


Track finding efficiency 


0.8 


0.6 


0.9 


0.3 


7r° finding efficiency 


0.0 


0.0 


0.0 


4.3 


Kg finding efficiency 


0.0 


0.0 


1.8 


0.0 


Unused tracks 


0.3 


0.3 


0.3 


0.3 


Signal shape fit function 


0.3 


0.3 


0.3 


4.2 


Simulation of FSR 


0.4 


0.4 


0.4 


0.4 


Simulation of form factors 


0.0 


0.1 


0.1 


0.3 


Limited MC statistics 


0.1 


0.2 


0.2 


0.3 


Total unc(M;lain(A- 


1.0 


1.1 


2.1 


6.1 



tector material (0.4%) is estimated by varying the amount of FSR modeled by the PHOTOS 
algorithm [37] and by repeating the analysis without recovery of photons radiated by the 
positron and comparing to the standard results. The uncertainty associated with the simu- 
lation of ISR (e+e" — > DD^y) is neghgible. There is a systematic uncertainty due to finite 
MC statistics (0.1% to 0.3% depending on mode). 

Table VI is a summary of the systematic uncertainties associated with the measurement 
of the four absolute semileptonic branching fractions. These estimates of systematic uncer- 
tainty are added in quadrature to obtain the total systematic uncertainty: 1.4%, 6.1%, 1.5%, 
and 2.4% for Tr-e+Ue, D+ 7r°e+i/e, ^° ^ K-e+Ue, and D+ K'^e+Ue, respectively. 

C. Comparison to Previous Measurements 

The branching fraction measurements with all tag modes combined for each of the 
four semileptonic modes reported in Table III, are in good agreement with previous 
CLEO-c measurements using the same technique [31] obtained with a smaller data sam- 
ple, and supersede them. In Table III we also compare our measurements to PDG 2004 [32] 
averages. We compare to PDG 2004 because subsequent PDG averages [38, 39] are domi- 
nated by our previous CLEO-c measurements. In Table VII we compare our measurements 
of B{D° — > K~e^Ue) and B{D^ — > n'e'^i'g) to previous measurements and to theoreti- 
cal predictions. Our measurements agree well with previous measurements including the 
CLEO-c neutrino reconstruction analysis [3], which we denote by "untagged" hereinafter. 

The widths of the isospin conjugate exclusive semileptonic decay modes of the 
and are related by isospin invariance of the hadronic current. The ratio r(D° — > 
K~e'^Ue)/r{D'^ — > K^e^Ue) is expected to be unity, while the corresponding ratio for pions 
is expected to be two. Using our results and the lifetimes tdo — 410.3(1.5) x 10~^^sand 
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TABLE VII: Comparison of B{D^ K~l^V() and B{D^ 7T~£^i'i) values among different ex- 
periments and theoretical predictions. The first uncertainty is statistical, the second is systematic. 
The third uncertainty in the BABAR measurement is from the normalization to B{D^ — > K~Tr'^). 





i<:-£+f£(%) 


TT-^+I/^ (0.1%) 


PDG (2004) [32] 


3.58(18) 


3.6(6) 


BES II (e) [40] 


3.82(40) (27) 


3.3(13)(3) 




O.t i \ i'±) 




LQCD (Abada) [24] 


2.99(45) 


2.4(6) 


QCD SR (Ball) [21] 


2.7(6) 


1.6(3) 


LCSR (KRWWY) [22] 


3.6(14) 


2.7(10) 


LCSR (WWZ) [41] 


3.9(1.2) 


3.0(9) 


CLEO-c (e) [31] 


3.44(10)(10) 


2.62(25)(8) 


Belle {e,n) [14] 


3.45(7) (20) 


2.55(19)(16) 


BABAR (e) [10] 


3.522(27) (45) (65) 




CLEO-c (tagged, e) 


3.61(5)(5) 


3.14(13)(4) 


CLEO-c (untagged, e) [3] 


3.56(3)(9) 


2.99(11)(9) 



= 1040(7) X lO-iSg [32]^ obtain 

and 



r(L>o 



= 1.03(2)(2) (22) 



1.04(9)(6), (23) 



2T{D+ TT^e+Ue) 

where correlated and uncorrelated systematic uncertainties are taken into account. These 
ratios are consistent with isospin predictions, and supersede the corresponding ratios in 
Ref. [31], which were measured with the same technique. These ratios are also consistent 
with the CLEO-c untagged analysis [3], and two less precise results: a measurement from 
BES II using the same technique [42] and an indirect measurement from FOCUS [43]. 

As the data are consistent with isospin invariance, the precision of each branching frac- 
tion can be improved by averaging the and results for isospin conjugate pairs. For 
the isospin-averaged semileptonic decay widths, with correlations among systematic uncer- 
tainties taken into account, we find 

r{D Ke+Ue) = 8.73(9) (15) x 10" V"^ (24) 

and 

r{D 7re+i/e) = 0.76(3)(2) x 10-V"\ (25) 

where for the latter partial width we have used r(£)° — > n'e'^Ue) — 2r(D+ — > Tr^e+i/g). 
The measured ratio of decay widths for D — > Tre+i/e and D — > Ke'^Ue provides a test of 
the LQCD charm semileptonic rate ratio prediction [26]. Using the results obtained in this 
analysis, we find 



0.0868(38) (4) (26) 
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TABLE VIII: Ratios of semileptonic decay widths of D and to the pseudoscalar mesons vr and 
K (first four hnes) and the isospin averaged ratio of semileptonic decay widths (fifth hne). The 
uncertainties are statistical and systematic. 



Ratios Measured values 

r{D^ ^ K-e+iye)/T{D+ K^e+iye) 1.030(24)(20) 
r(D° ^ 7r-e+i/e)/2r(D+ ^ 7r°e+i/e) 1.037(86)(57) 
r(£)° ^ 7r-e+i/e)/r(D° ^ -fi:~e+i/e) 0.0868(38)(4) 
2r(D+ ^ 7rOe+i/e)/r(D+ ^ ^Oe+z/g) 0.0863(64)(53) 
T{D 7re+iye)/T{D Ke+u^) 0.0868(33) (14) 



and 

H^^^^^i^^ 0.0863(64,(53). (27) 

These results are consistent with LQCD [26] and with previous measurements [3, 12, 44]. 
Finally, by averaging the I?° and results for isospin conjugate pairs we obtain 

i^§^^^0.0868(33)(14), (28) 

where we have again used V{D'^ — >• 7r~e"'"i/e) = 2r(£)+ — >• Tr^e+i/e). A complete set of ratios 
of partial semileptonic decay widths measured in this analysis is given in Table VIII. 



VI. STUDY OF SEMILEPTONIC DIFFERENTIAL DECAY RATES 
A. Measurement of the Differential Decay Rate 

We now describe how the efficiency-corrected absolutely-normalized differential decay 
rate distributions are obtained. Full event reconstruction allows a direct measurement of 
the neutrino momentum with excellent resolution. The invariant mass squared of the e^i/e 
pair, g^, is calculated in the -0(3770) rest frame in the following way (using as an example 

= (^beam - Ek)^ - (-ptag - PK^, (29) 
Ptag =Ptag-\/^beam-"^Z), (30) 

where and px are the energy and three-momentum of the kaon. The q'^ resolutions 
(^reconstructed " ?|enerated) avcragcd ovcr the _entire g2 range are about 0.012 (GeV/c^)^ for 
TT'e+Ue, D° K'e+u^ and D+ K^e+u^, and approximately 0.040 {GeV/c^f for 
D+ 7r°e+i/e- For D+ n^e+Ue, the ^reconstructed - ^generated distribution is well described by 
a Gaussian. For other semileptonic modes these distributions are consistent with a double 
Gaussian with a's that differ by a factor of 2.5, with the wider Gaussian mostly due to FSR. 

As the D mesons are produced almost at rest at the ■ip{3770), and the CLEO-c detector is 
nearly hermetic, the semileptonic reconstruction efficiencies are almost constant across the 
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range. In consequence the shape of the spectrum receives only minor distortions due 
to detector acceptance. The excellent q^ resolution likewise leads to only minor distortions 
due to q^ smearing. 

Events satisfying the reconstruction criteria of Sec. IV that lie in the U signal region, 
defined as —60 MeV < f/ < 60 McV, are sorted into bins of q^. Ten bins of equal 
size (o'max/lO) used for K^e^Ug and K^e'^Ug. Nine (seven) bins are 

used for — > 7r~e"'"i/e {D~^ vr^e+t'e) with the last bin two (four) times wider than the 
other bins to allow for the smaller number of events at large q"^ for these modes. The bin 
limits are given in Table IX. 

TABLE IX: The upper edge of each bin in units of GeV^/c^ for each semileptonic mode studied 
in this work. 



Mode 


Bin 1 Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 Bin 7 Bin 8 Bin 9 Bin 10 




TT e'^l/e 


0.30 0.60 0.89 1.19 1.49 1.79 2.08 2.38 ql,^ 
0.30 0.60 0.90 1.20 1.50 1.80 ql,^ 
0.19 0.38 0.56 0.75 0.94 1.13 1.32 1.50 1.69 ql,^ 
0.19 0.38 0.56 0.75 0.94 1.13 1.32 1.51 1.69 ql,^ 



TABLE X: Numbers of events, estimated backgrounds and yields in q^ bins for the four semileptonic 
modes. The uncertainty in parentheses is statistical. The q^ bins are defined in Table IX. 



Mode 




Bin 1 


Bin 2 


Bin 3 


Bin 4 Bin 5 Bin 6 


Bin 7 


Bin 8 


Bin 9 


Bin 10 






Number of events 


130(11) 


122(11) 


99(10) 


105(10) 76(9) 56(8) 


66(8) 


38(6) 


19(4) 




DO-^ 


TT e+fe 


Background 


8.9(7) 


8.3(7) 


7.0(6) 


6.2(5) 4.7(4) 3.5(3) 


2.8(2) 


2.4(2) 


2.9(2) 








Yield 


121(11) 


114(11) 


92(10) 


99(10) 71(9) 52(8) 


63(8) 


36(6) 


16(4) 








Number of events 


48(7) 


46(7) 


44(7) 


36(6) 34(6) 30(6) 


48(7) 








D+ ^ 




Background 


1.8(1) 


1.6(1) 


2.5(2) 


3.0(2) 2.7(1) 3.1(2) 


20.0(1.4) 












Yield 


46(7) 


44(7) 


42(7) 


33(6) 31(6) 27(6) 


28(7) 












Number of events 


1239(35) 


1169(34) 1006(31) 923(30) 821(29) 594(24) 


464(22) 


293(17) 139(12) 


29(5) 




K-e+Ue 


Background 


6.7(6) 


6.7(6) 


8.1(7) 


7.7(7) 9.1(8) 8.7(7) 


5.3(5) 


3.9(3) 


3.1(3) 


1.5(1) 






Yield 


1232(35) 1162(34) 


998(32) 


915(30) 811(29) 585(24) 


459(22) 


290(17) 136(12) 


28(5) 






Number of events 


570(24) 


502(22) 


442(21) 


379(19) 298(17) 255(16) 


210(14) 


112(11) 


64(8) 


19(4) 


D+ ^ 




Background 


2.4(3) 


3.0(4) 


3.4(5) 


3.8(5) 3.3(4) 3.5(5) 


2.9(4) 


2.1(4) 


1.8(2) 


1.2(2) 






Yield 


568(24) 


499(22) 


439(21) 


375(19) 295(17) 251(16) 


207(15) 


110(11) 


62(8) 


17(4) 



The number of events in the data, the estimated background, and the background- 
subtracted yield in each bin of are provided in Table X. To obtain dV /dq^ for each 
semileptonic mode, the background is subtracted from the observed distribution. The 
number of signal events N^^^'^^ in the ith bin is given by 

^tag,SL _ ^ ^ij ^produced ^g-j^'^ 
3 

where e*-* is the semileptonic efficiency matrix which accounts for acceptance and resolution 
effects. This matrix equation is inverted to obtain ]\jY°'^^^'^^ a vector of efficiency corrected 
signal events with a D tag in the data. When properly normalized, the elements of ]S[Y°'^^^'^ 
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give the absolute decay rate in bins. Efficiency matrices, e*^, for eacli semileptonic mode 
are obtained using signal MC samples. The procedure for calculating the efficiency matrices 
is analogous to that for e: 

tig,SL/^tag, (32) 



^3 — ^0 



with e^-?^g SL obtained as 

M _ signal ioo\ 
^tag.SL - > W-^J 

where A^^tai ^'^^ number of signal events generated in the jth bin, and A^g/gnai ^^^^ 
number of signal events that are generated in the jth bin and reconstructed in the ith 
bin. Efficiency matrices for each of the four modes are given in Table XI. These efficiency 
matrices have been calculated for the simple pole model, with the distribution re- weighted 
to the value of Mpoic determined by the data for each mode, and weighted by the tag yields 
given in Table I and II. We note that at the present level of precision, due to the use of 
efficiency matrices combined with the fine binning in g^, the values we determine for the 
shape and normalization parameters in the form factor fits are not sensitive to the model used 
to generate the efficiency matrices. The statistical uncertainty of the background-subtracted 
and efficiency-corrected decay rate distribution for each bin is given by 



2 



+ [^(e-^)]?[iV4g,SL]')- (34) 



TABLE XI: Abridged efficiency matrices for the four semileptonic modes. Matrix elements are in 
percent. The semileptonic efficiency inchidcs subsidiary branching fractions [36]. Those that are 



not on or adjacent to a diagonal are null and are not shown. The uncertainties are statistical. The 
bins are defined in Table IX. 





e{i,j) e{l,j) e{2,j) e(3,j) e{4,j) e{5,j) e{6,j) e{7,j) e{8,j) e{9,j) e{lQ,j) 




e{i,i-l) 0.0(0) 2.08(5) 1.98(5) 1.85(5) 1.67(5) 1.53(6) 1.23(6) 1.21(6) 0.85(7) 
e(i,i) 62.81(26) 64.13(28) 67.25(31) 69.10(,34) 69.48(37) 70.36(42) 70.23(49) 70.02(61) 68.17(75) 
e{i,i + l) 2.23(5) 2.14(5) 2.00(6) 1.72(6) 1.58(6) 1.45(7) 1.21(8) 0.93(7) 0.0(0) 


D+ -» 7r°e+i/e 


e{i,i-l) 0.0(0) 3.32(9) 3.68(10) 3.35(10) 2.82(10) 2.62(11) 2.50(12) 
e{i,i) 36.00(29) 35.48(30) 36.17(33) 36.55(36) 35.84(39) 35.67(44) 38.42(37) 

e{i,i + l) 2.21(7) 2.07(8) 2.06(9) 1.87(9) 1.72(10) 0.81(5) 0.0(0) 


D° -* K-e+Ue 


e{i,i-l) 0.0(0) 2.40(3) 2.36(4) 2.24(4) 1.99(4) 1.81(4) 1.48(4) 1.16(4) 0.85(5) 0.47(2) 

52.86(16) 52.78(17) 55.56(18) 57.98(20) 59.49(23) 59.72(25) 58.96(30) 57.43(36) 53.78(50) 39.62(86) 
e{i,i + l) 2.70(4) 2.60(4) 2.54(4) 2.37(5) 2.33(5) 2.07(6) 1.86(7) 1.63(9) 1.52(17) 0.0(0) 


D+ K"( +Ur 


e(i,i-l) 0.0(0) 0.82(2) 0.75(2) 0.73(2) 0.64(2) 0.62(2) 0.55(2) 0.49(2) 0.38(2) 0.21(2) 

<-(?.?■) 18.14(8) 17.57(8) 18.12(8) 18.64(9) 18.63(10) 20.20(11) 19.01(13) 19.17(17) 20.12(24) 20.13(51) 
(((./ + 1) 0.81(2) 0.82(2) 0.82(2) 0.79(2) 0.72(2) 0.03(2) 0.(i()(3) 0.03(1) 0.78(10) 0.0(0) 



Background-subtracted, efficiency-corrected and absolutely normalized decay rate distri- 
butions for the four semileptonic modes are given in Table XII. They constitute the main 
result of this analysis and can be used to compare to other experimental measurements and 
to theory without a need for knowledge of CLEO-c acceptance and resolution. 

Table XII includes statistical uncertainties and the associated correlation matrices. As 
discussed in Sec. VI C 1, systematic uncertainties are approximately fully correlated between 

bins across the entire q^ range. Therefore we include systematic uncertainties for each 
bin for each semileptonic mode in Table XII without correlation matrices. 
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TABLE XII: Absolutely normalized decay rates in bins of q^, with statistical and systematic 
uncertainties in parentheses, are given in the first row for each decay mode. These distributions 
are background-subtracted and efficiency-corrected. The truncated statistical correlation matrices 
are shown in the last three or five rows. The bins are defined in Table IX. 



Bin: 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


r(7r-e+i/e) 


1.482 


1.325 


1.014 


1.087 


0.774 


0.560 


0.697 


0.388 


0.183 




[ns-i] 


(149)(24) (142)(24) (122)(19) (121)(21) (103)(15) 


(87)(11) 


(94)(13) 


(72) (7) 


(55) (4) 








-0.068 


-0.063 


-0.056 


-0.049 


-0.045 


-0.039 


-0.036 


-0.027 






1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 






-0.068 


-0.063 


-0.056 


-0.049 


-0.045 


-0.039 


-0.036 


-0.027 








1.432 


1.270 


1.167 


0.908 


0.917 


0.799 


0.786 








[ns-l] 


(233) (97) (233) (77) (226) (63) (204) (44) (202) (39) (192) (30) (253) (25) 








Cu—2 






0.020 


0.019 


0.015 


0.011 


0.002 








'^11— I 




-0.155 


-0.162 


-0.154 


-0.130 


-0.123 


-0.078 








Cii 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 








Cii+i 


-0.155 


-0.162 


-0.154 


-0.130 


-0.123 


-0.078 










Cii+2 


0.020 


0.019 


0.015 


0.011 


0.002 












r{K-e+Ue) 


1.767 


1.601 


1.303 


1.160 


1.015 


0.725 


0.586 


0.383 


0.193 


0.053 


[10 X ns-i] 


(53) (27) 


(52) (26) 


(46) (22) 


(42) (20) 


(39) (18) 


(33) (13) 


(29)(11) 


(24)(7) (18)(4) (11)(1) 


Cii— 1 




-0.096 


-0.092 


-0.084 


-0.074 


-0.069 


-0.059 


-0.051 


-0.043 


-0.036 


Cii 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


Cii+1 


-0.096 


-0.092 


-0.084 


-0.074 


-0.069 


-0.059 


-0.051 


-0.043 


-0.036 




r(/?Oe+i.e) 


1.785 


1.540 


1.320 


1.105 


0.870 


0.757 


0.613 


0.312 


0.175 


0.050 


[10 X ns-l] 


(79) (42) 


(77) (38) 


(70) (33) 


(63) (29) 


(55)(23) 


(52)(20) 


(46)(17) 


(33)(9) 


(24) (5) 


(13)(2) 


Cii— 1 




-0.093 


-0.088 


-0.084 


-0.077 


-0.072 


-0.063 


-0.060 


-0.049 


-0.040 


Cii 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


Cii+1 


-0.093 


-0.088 


-0.084 


-0.077 


-0.072 


-0.063 


-0.060 


-0.049 


-0.040 





TABLE XIII: Form factor distributions in bins of in the first row, with statistical and systematic 
uncertainties in the second row in parentheses for each decay mode. The bins are defined in 
Table IX. The entries for ir^e'^Ue have been scaled by the isospin factor \/2. 



Bin: 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


fl{q^)\Vcd\(n-e+Ue) 


0.160 


0.175 


0.180 


0.222 


0.230 


0.249 


0.370 


0.398 


0.458 






(8)(1) 


{9)(2) 


(11)(2) 


(12)(2) 


(15)(2) 


(19)(2) 


(25)(4) 


(37) (38) 


(68)(5) 




fl{q'')\Vcd\{7T°e+Ue) 


0.156 


0.170 


0.191 


0.202 


0.249 


0.296 


0.357 










(12)(5) 


(15)(5) 


(18)(5) 


(21)(5) 


(26)(5) 


(33) (6) 


(54) (6) 








f^(q^)\Vcs\(K-e+u,) 


0.759 


0.806 


0.821 


0.887 


0.968 


0.980 


1.098 


1.180 


1.268 


1.519 




(12)(6) 


(14)(6) 


(15)(7) 


(17)(8) 


(20)(9) 


(24) (9) 


(30)(10) 


(40)(11) 


(63) (13) 


(159)(19) 


f^(q')\Vr,\(K'\'+,^^) 


0.760 


0.788 


0.821 


0.862 


0.893 


0.997 


1.118 


1.061 


1.200 


1.437 




(18)(<)) 




[2:i)(W) 


[mill) 


m(i2) 


(:-!(,) (l:-!) 


(11) (15) 


(.-,8)(15) 


(87)(1V) 


(202) (25) 



The partial differential decay rates are expected to be identical by isospin invariance. 
A powerful check of our understanding of the data is therefore provided by comparing the 
background-subtracted, efficiency-corrected rates in Table XII for D'^ and D^. We make 
the comparison by removing the kinematic term and constants from the differential rate to 
reveal /+(5^)|V^|, where q — d or s, 



eg I 



dV 24773 



(35) 



\dq^Glp%,' 

where dT/dq^ is obtained by dividing the integrated rate in each q^ interval by the corre- 
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spending bin size, and in the ith q bin is given by 



I f -\- (^center of bin i 



(36) 



where the form factor parameters are measured in the data using the three parameter series 
parametrization (see Sec. VIB). For D Ke^v^ {D Tre+z/g), /+(q'^)|K?| varies by only 
a factor two (three) across the range. Table XIII and Fig. 4 show /+(g^)|V^| and f+{q^) 
in data for all four semileptonic modes. The isospin conjugate distributions are consistent. 



2.5 
2.0 
1.5 
1.0 



o D° ^ 71" e+ 



+0 



0120708-004 



-I I I I I I I I l_ 







1.0 



2.0 0.5 

q2 (GeV2/c4) 



1.0 



1.5 



1.5 



1.0 



FIG. 4: The data displayed as /+((?^) for the four semileptonic modes. In each case f+{q^) is 
absolutely normalized. The D'^ and distributions are each offset symmetrically in to facilitate 
display. 



B. Fitting the Differential Decay Rate to Determine Form Factors 

We use the least squares method to fit the absolutely-normalized efficiency-corrected and 
background-subtracted distributions. A is constructed from differences between the 
number of efficiency corrected signal events with a D tag in the ith (f bin, -/Vproduced> ^^"^ 
the theoretically predicted number of events in the ith q^ bin, -/Vpredicted' ^ given set of 
form factor parameters, where -^Vpredicted obtained using 



predicted 



NugTD \f+[q ,e)\ dq , (37) 



and where td is the lifetime of the relevant D meson, and 6 is the vector of form factor 
parameters that govern the decay rate. Taking into account the correlations among the bins 
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and the correlations among the elements of the inverted efficiency matrix [45], the is 
given by 




(38) 



where A,-, is 



^ij — X 



tagjSL 



(39) 



k 



Systematic uncertainties and correlations among them are not included in the fit. Instead 
a systematic uncertainty from each source is estimated separately The fitting procedure has 
been tested using ensembles of fits to 100 mock data samples that each correspond to the 
same integrated luminosity as the data, for a wide range of values of the form factor param- 
eters. It has been established that the statistical Tinccrtaintics from this fitting procedure 
are consistent with the smallest statistical uncertainties expected from a fit, estimated using 
the Cramer - Rao inequality [46], and that the fit is consistent with being unbiased. 

Fits to the data are made for two parameters related to the shape and the normalization 
of the f+{q'^) form factors for the series parametrization, the simple pole, the modified pole 
model, and the 1SGW2 model. For the series parametrization we also present results of 
fits for three parameters, where the third parameter is a second shape parameter. As an 
example. Fig. 5 shows simultaneous fits to modes related by isospin symmetry. Before pre- 
senting numerical results of the form factor measurements, we describe a study of systematic 
uncertainties in the next section. 

C. Study of Systematic Uncertainties for dT/dq^ and Form Factor Measurements 

1. Systematic Uncertainties for dT / dq^ in bins 

Each source contributing systematic uncertainty to the absolute branching fractions also 
contributes systematic uncertainty to measurements of the partial rate dV /dq^ in bins. 
Procedures identical to those used in the absolute branching fraction measurements are em- 
ployed to estimate systematic uncertainties for dT/d(f. In addition, there is a systematic 
uncertainty associated with imperfect knowledge of D'^ and meson lifetimes. Table XIV 
reports total systematic uncertainties, and the separated correlated and uncorrelated com- 
ponents, for dV /dq^ in q^ bins for the four scmilcptonic modes. 

Systematic uncertainties associated with finding and identifying the hadron (positron) in 
the final state of a semileptonic decay are measured in bins of hadron (positron) momentum 
and propagated to the dV /dq^ distributions. In the rest frame of the decaying D meson, 
is determined by the momentum of the final state hadron. Because ■?/'(3770) decays produce 
D mesons with a small boost, q^ is strongly correlated with the momentum of the final state 
hadron measured in the laboratory frame. Therefore systematic uncertainties measured in 
hadron momentum bins, when propagated to dT/dq^, lead to uncertainties that are mostly 
uncorrelated between bins. The correlation between and the positron momentum in 
the laboratory frame is less pronounced due to additional degrees of freedom associated with 
the undetected neutrino. Systematic uncertainties in positron momentum bins are therefore 
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0120708-005 




q2 (GeV2/c4) 

FIG. 5: Simultaneous three parameter fits to the background-subtracted, efficiency-corrected and 
absolutely-normalized data derived decay rates for modes related by isospin D ■K~e~^i'e; 
7r°e+fe, and — K~e^Ue\D^ — K'^e^Ue- The — 7r°e+i/e distribution is scaled by a factor 
two to account for isospin. The series parametrization with a quadratic term is used for these fits. 
The fit result is the line. 

averaged over a range in and their net effect is to produce uncertainties that are nearly 
constant and fully correlated between bins. 

To simplify the estimation of the systematic uncertainties, we assume that a given sys- 
tematic uncertainty is either fully correlated or uncorrelated between (f bins as discussed in 
the remainder of this section. 

Studies of the momentum dependence of the systematic uncertainty associated with track 
finding efficiencies are performed in three momentum bins covering the entire momentum 
range accessible in D meson decays at the -0(3770). Efficiencies for positively and negatively 
charged pions and kaons are measured separately. We assume that track finding efficiencies 
for positrons are identical to those for positively charged pions. A systematic uncertainty 
from track finding efficiencies in a bin is calculated by weighting charged hadron (positron) 
spectra with the efficiency uncertainties measured in the hadron (positron) momentum bins 
and summing contributions from different momentum bins. Due to the coarse binning used 
in the tracking studies and because positron and charged hadron track finding uncertainties 
are combined in each q^ bin, systematic uncertainties associated with track finding efficiency 
are strongly correlated between q^ bins. We assume that they are fully correlated. 

Systematic uncertainties associated with charged hadron identification in K'e^i/f. 
{D^ Ti~e^V(.) are obtained by weighting charged hadron spectra with the statistical uncer- 
tainties associated with charged hadron identification measured in 100 (80) MeV/c - wide 
momentum bins and summing contributions from different momentum bins in quadrature. 
Because the hadron momentum is strongly correlated with (f, these systematic uncertain- 
ties are largely independent for well-separated values of q^. We therefore assume that the 
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TABLE XIV: Correlated, uncorrelated and total systematic uncertainties in dV/dq'^ for each bin 
for the four semileptonic modes. The q^ bins are defined in Table IX. 



Mode 



Type 



Uncertainty (%) 
Bin 1 Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 Bin 7 Bin 8 Bin 9 Bin 10 







Correlated 


1 

i.O 


i. i 


1 c 
i.O 


1 o 

i.y 


1 Q 
i.O 


1 c 
i.O 


i.O 


1 c 
i.O 


z.U 




JJ IT 




Uncorrelated 


U.O 


U.4 


0.5 


0.5 


O.D 


O.D 


O.D 


0.8 


1.1 








Total 


1.6 


1.8 


1.9 


1.9 


1.9 


1.9 


1.9 


1.9 


2.3 








Correlated 


6.8 


6.1 


5.4 


4.8 


4.2 


3.6 


3.1 








D+ ^ n^e+ 




Uncorrelated 


0.7 


0.8 


0.8 


0.9 


0.9 


1.0 


0.9 












Total 


6.8 


6.1 


5.4 


4.9 


4.3 


3.8 


3.2 












Correlated 


1.4 


1.5 


1.6 


1.6 


1.7 


1.8 


1.8 


1.7 


1.9 


1.5 


K-e 




Uncorrelated 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.5 


0.7 


0.9 


1.9 






Total 


1.5 


1.6 


1.7 


1.7 


1.8 


1.8 


1.9 


1.9 


2.1 


2.4 






Correlated 


2.4 


2.5 


2.5 


2.5 


2.6 


2.6 


2.6 


2.7 


2.6 


2.6 


D+ K^e 




Uncorrelated 


0.4 


0.4 


0.5 


0.5 


0.6 


0.6 


0.7 


0.9 


1.2 


2.2 






Total 


2.4 


2.5 


2.5 


2.6 


2.6 


2.7 


2.7 


2.9 


2.9 


3.4 



systematic uncertainties associated with hadron identification are uncorrelated between (f 
bins. 

The systematic uncertainty in the 7r° reconstruction efficiency varies from 1.3% for low 7r° 
momenta to 6.3% for high 7r° momenta and is found to be fully correlated between bins. 
Systematic uncertainties associated with the reconstruction are found to be independent 
of the momentum and are fully correlated between bins. 

Systematic uncertainties due to simulation of ISR and FSR are strongly correlated be- 
tween (f bins. Systematic uncertainties from FSR are assigned based on differences between 
the main results and results of fits with efficiency matrices obtained using a subset of signal 
MC events without FSR. To evaluate systematic uncertainties associated with ISR, we re- 
peated the analysis with two alternative efficiency matrices: one using signal MC events with 
soft ISR photons {E^ < 25 keV) and the other from the remainder of the signal MC events. 
Comparing results of fits with these two efficiency matrices, we conclude that systematic 
uncertainties due to ISR are negligible. 

The background is modeled using the generic MC sample. Systematic uncertainties as- 
sociated with the modeling of background are obtained by varying the composition of the 
background sample according to uncertainties in the branching fractions of processes produc- 
ing background, and by the statistical uncertainties in the normalization for each background 
component. In addition, in cases where a background component arises from misidentified 
hadrons or Icptons, background normalizations are varied according to the uncertainty in 
the relative misidcntification rates between the data and MC simulation. 

Systematic uncertainties from imperfect knowledge of the (0.4%) and (0.7%) 
meson lifetimes, the number of tags (0.4%), and unused tracks (0.3%) are fully correlated 
between bins. Systematic uncertainties due to the limited size of the MC samples used 
to measure the efficiency matrices are statistical in origin and are therefore uncorrelated 
between q^ bins. 
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Three systematic uncertainties for each bin are presented in Table XIV. These are the 
combined sum in quadrature of all correlated and all uncorrelated contributions, and the 
total systematic uncertainty. The magnitude of the systematic uncertainty in each bin is 
significantly smaller than the corresponding statistical uncertainty, and the relative size of 
the uncorrelated systematic uncertainty is small compared to the correlated systematic un- 
certainty in nearly all bins. (Note, in the last q^ bin the uncorrelated systematic uncertainty 
is dominated by uncertainty due to the limited size of the MC sample, and for D — > Ke^u^ 
it is comparable to the correlated systematic uncertainty.) For comparison to theory and 
for the form factor measurements presented here we assume that systematic uncertainties 
are fully correlated between bins. 

2. Systematic Uncertainties for Measurements of f+{0) and Form Factor Shape Parameters 

The normalization parameter, /_|_(0)|V^g|, and form factor shape parameters are deter- 
mined from simultaneous two parameter fits to dV /dq^{(f) for each isospin conjugate semilep- 
tonic mode. In each case the correlation coefficient between the form factor shape parameter 
and the normalization parameter is found to be small. 

Systematic uncertainties associated with the absolute form factor normalization, /+(0), 
for each semileptonic mode, are one half the systematic uncertainties in the branching frac- 
tion measurements presented in Sec. V B combined in quadrature with the small uncertain- 
ties associated with the knowledge of (0.4%) and D+ (0.7%) [32] lifetimes and the CKM 
matrix elements \Vcs\ (0.1%) and \Vcd\ (1-3%) obtained from the unitarity constraints of the 
CKM matrix. 

Systematic uncertainties for form factor shape parameters are obtained from one param- 
eter fits with absolute normalizations fixed. In the rest of this section, we describe sources 
of systematic uncertainty for form factor shape parameters and how they are estimated. 

A systematic uncertainty associated with the fit procedure is assigned by examining 
the pull distributions resulting from fits to ensembles of mock data samples. The studies 
determine that the fit has good fidelity, and place an upper limit on the existence of bias at 
15% of the statistical uncertainty in the measurement on data, which we take as a systematic 
uncertainty associated with the fit method (Table XV). 

As discussed in the previous section, most sources of systematic uncertainty are (f inde- 
pendent and, consequently, do not contribute a systematic uncertainty in the form factor 
shape parameter measurement. Accordingly, to assign systematic uncertainties for tracking 
efficiency, and 7r° finding, and hadron and electron identification, a correlation with par- 
ticle momentum consistent with our knowledge of each systematic effect is introduced. This 
is achieved by constructing a model according to which each systematic uncertainty varies 
linearly as a function of the particle momentum. The slope for each systematic uncertainty 
is determined by the precision with which each systematic effect is known. We fit the data 
using efficiency matrices modified according to this model. We also construct a set of mock 
data samples for each source using the model and fit them to obtain systematic uncertainties 
for the form factor shape parameters. We find that systematic uncertainties measured by 
these two methods are consistent. 

Systematic uncertainties associated with the simulation of FSR and from background es- 
timation are obtained as described in the previous section. The total systematic uncertainty, 
given in Table XV, ranges from 19% to 53% of the statistical uncertainty. The ratio of the 
systematic to statistical uncertainties for shape parameters are found to be consistent for 



27 



all parametrizations. 

TABLE XV: Systematic uncertainties for form factor shape parameters in units of the statistical 
uncertainty of the measurement with data. The modes are labeled by their final state hadrons. 



Systematic uncertainty {(Jstat) 



Sources 


K- 




TT 


ttO 


Track finding 


0.03 


0.00 


0.02 


0.00 


Kg finding 


0.00 


0.20 


0.00 


0.00 


ir^ finding 


0.00 


0.00 


0.00 


0.38 


Hadron ID 


0.11 


0.00 


0.05 


0.00 


Electron ID 


0.02 


0.01 


0.01 


0.01 


FSR 


0.20 


0.16 


0.07 


0.04 


Background 


0.05 


0.07 


0.06 


0.34 


MC size 


0.10 


0.10 


0.05 


0.05 


Fitter 


0.15 


0.15 


0.15 


0.15 


Total 


0.30 


0.32 


0.19 


0.53 



D. Form Factor Measurement Results 

The fit described in Sec. VI B is appfied to the decay rates in Table XII for each 
semileptonic mode and to pairs of modes related by isospin. Five fits are carried out per 

mode. In each case the normalization parameter /+(0)|V^q| and one or more form fac- 
tor shape parameters are determined. Specifically the shape parameters are Mpoie (simple 
pole model), a (modified pole model), and r (1SGW2). For the series parametrization 
we map the data to the variable z. The quantity P{z)(p{z)f^{z) is, by convention, con- 
strained to unity ai z = z^^^x, which corresponds to = 0. We fit to the distribution 
dr/dq^ = {Gl/247r^)p^j^jVcg\^alFl{z), where: F+{z) = [P{z)(l){z)]-^ x (l + n^ + ra^^) with 
ri = ai/oo, r2 = Oa/oo, and /+(g^) = aoF+{z). 

The fit returns the normalization parameter /_|_(0)|Vcg| and either ri or ri and r2- We 
test the sensitivity of the data to the number of parameters and the convergence of the 
series. For the series parametrization the slope at the intercept, 1 + 1/(3 — 6, is also reported. 
Results of fits to each parametrization are given in Table XVI and Table XVll. 

Comparisons of four of the five fits to the data for each of the four modes are shown in 
Fig. 6 (ISGW2 is excluded). To facihtate a comparison, in Fig. 7 we normalize each fit to the 
result of the three parameter series fit. It can be seen that each of these parametrizations 
provides an adequate, and almost identical, description of the data when the shape parameter 
is allowed to be free. To illustrate the difference between the linear and quadratic ;z-expansion 
fits. Fig. 8 shows P{z)4>{z)f+{z)/P{z{q'^ = 0))((){z{q'^ = 0)) for both as a function of z. 

An independent assessment of the quality of the fits to the data is obtained from the 
ability of the fit to describe distributions in the data in two variables that are not used to 
constrain the fit. The first variable is the angle between the in the D meson frame and 
the positron in the frame, Owe- The second variable is the laboratory momentum of the 
positron, \pe\. Figure 9 shows distributions for cos 9we and \pe\ in data and the projections 
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FIG. 6: Projections of fits onto f+{q^) for each semileptonic mode. In each case f+{q^) is absolutely 
normalized. The data are shown as points with error bars. The lines are fits to the simple pole 

model (long dash), the modified pole model (short dash), the series parametrization with two free 
parameters (dot), and the series parametrization with three free parameters (solid). 



of the fit, where the background contributions are shown as hatched histograms. The fits 
describe the distributions in these two variables well. 

Using the ISGW2 parametrization we determine the isospin conjugate average values of 
the meson radius to be 

= 1.53(4)(2) GeV-\ (40) 

r'^ = 1.95(10) (5) GeV-^ (41) 

They are the most precise measurements of these quantities to date, and are 18cr and 5a 
from the ISGW2 expected values = 1.12 GcV^^ and r'^ = 1.410 GcV~^, respectively. 
We have assigned no uncertainty to the theoretical prediction, and assume here and in what 
follows, that the experimental uncertainties derived from the fit are Gaussian distributed. 
A comparison to other recent measurements is given in Table XVIII. The measurements by 
BABAR and Belle disagree by 3.8(7. Our measurement of is over 4(7 smaller than Belle, 
and 1.6(T smaller than BABAR. 

Using the simple pole model, we determine the isospin conjugate average pole masses to 

be 

Mp^i, = 1.97(3)(1) GeV/c^ (42) 
M;„i, = 1.95(4)(2)GeV/c2. (43) 
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FIG. 7: Form factor fit comparison for each semileptonic mode. All data (squares) and fits (his- 
tograms) arc normalized to the relevant three parameter fit result (Series (3) solid line at unity). 
The simple pole, modified pole, and two parameter series fit (Series (2)) are represented by long 
dash, short dash, and dotted histograms, respectively. 



These values differ by 4.6(7 and 1.3a from the well-measured masses: M^*+ = 2112.0 ± 
0.6 MeV/c^ and M^d.+ = 2010.0 ± 0.4 MeV/c^ [38], respectively. Comparison to previous 
measurements arc given in Table XIX and Table XX. For D — > Ki^u^ all of the more recent 
measurements are below the mass of the D* meson. Our measurement of M^^i^ is in excellent 
agreement with Ref. [3], but is 2.3a larger than the BABAR [10] measurement, and 2.5a 
larger than the Belle [14] measurement. For D — > ni'^u^ all measurements are much less 
precise and are in reasonable agreement, albeit within large uncertainties. All measurements 
are below the mass of the D* meson. 

Using the modified pole model, we determine the isospin conjugate average shape param- 
eters to be 

- 0.21(5)(2), (44) 

0.16(10)(5). (45) 

The values of and are 27a and 11a, respectively, from the values of ~ 1.75 and 
1.34 required by the BK parametrization. A comparison to previous measurements is 
given in Tabic XXI. For Ke'^Ug there is excellent agreement between this result and 
CLEO-c (untagged) [3], good agreement with previous measurements by CLEO III [12], and 
FOCUS [13], and QCD sum rules [21], but our result is lower than BABAR [10] by 2.6a, 
lower than Belle [14] by 2.7a and lower than the LQCD fit [26] by 4.2(7. The significance 
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FIG. 8: Series paramctrization form factor fit comparison to the data for two and three expansion 
parameters for each semileptonic mode. Series(2) (black line) with ztlo" uncertainty (black dashed 
lines). Series (3) (gray line) with ztla uncertainty (gray dashed lines). 



of the discrepancy between our result and the LQCD fit cannot be quantified rigorously, as 
the covariance matrix for the LQCD form factor is lost during the chiral extrapolation [26] . 
For Tre'^i'e there is reasonable agreement with CLEO-c (untagged) [3] and other previous 
measurements, albeit within large uncertainties. Our measurement of a'^ is 2.4cr smaller 
than the LQCD fit. 

Fits to the data using the first two terms of the z expansion are reported in Table XVL 
Fits using the first three terms are given in Table XVII and shown in Fig. 8. The expansion 
parameters are not predicted. The central value of the ratio of expansion parameters r2 is an 
order of magnitude larger than ri, however the statistical uncertainty is of similar magnitude 
to the central value, and therefore no statement can be made about the convergence of the 
expansion. Moreover, the data lack the precision, even in the copious D — > Ke'^Ue mode, 
to determine r2. For this reason there is no appreciable difference between the probability 
of the between the two parameter series expansion and three parameter series expansion 
fits for any mode. The compatibility of the data with linear dependence is consistent with 
the modified pole ansatz for f+{q^). Recently BABAR [10] using a data sample of 75,000 
— > K'e^fe events, found r\ — —2.5(2) (2) and — 0.6(6) (5), and that the differential 
rate is well-described by the z expansion with only a linear term. The results reported 
here for ri and r2 are in excellent agreement with CLEO-c (untagged) [3] and agree with 
BABAR [10] to better than 2cr with the precise level depending on the correlation coefficient 
for the BABAR ri and r2 parameters. 

The quadratic series expansion fit returns isospin conjugate average values for l-\-l/(5 — 5 
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FIG. 9: Distributions of positron momenta (left) and cos f?p^'e (right) in data (points with error bars) 
and the projection of the fit (sohd hne) for each semileptonic mode. The background contributions 
are the hatched histograms. 



of 0.69(12) (5) and 0.98(23) (10) for D Ke+v^ and D ne+iy^, respectively. These values 
are lOcr and 4a from the value of ~ 2 required by the BK parametrization, and are consistent 
with the results in [3] given in Tabic XXll. 

When the shape parameters are not fixed the parametrizations of the simple pole 
model, the modified pole model, the ISGW2 model, and the series expansion with two and 
three parameters are functionally almost identical over the range accessible in D meson 
semileptonic decay. For this reason each parametrization is able to describe the data with a 
comparable probability. 

Measurements of /+(0)|V"cg| are given in Table XVI for the 1SGW2, simple pole, modified 
pole, and two parameter series parametrization, and in Table XVII for the three parameter 
series parametrization. As each parametrization is able to describe the data, measure- 
ments of /+(0)|V^g| are very similar among parametrizations. For D Ke^Ue the values of 
/:f"(0)|V'cs| span about one half of a statistical sigma between the pole model, modified pole 
model, and series expansion (linear). However, the fit to the series expansion including a 
quadratic term returns a value of /:f^(0)|T4s| one statistical sigma larger than for the series 
expansion using a linear term. The statistical uncertainty is also increased by one third. 

For D — > ne'^i'e the values of /^(0)|V'cd| span a statistical sigma among the pole model, 
modified pole model, and the series expansion (linear). The fit to the series expansion in- 
cluding a quadratic term returns a value of (0)114^1 that only differs in the least significant 
digit from the value obtained for the series expansion using a linear term, but the statistical 
uncertainty is increased by one third. 

Using \Vcs\= 0.97334±0.00023 and \Vcd\ = 0.2256±0.0010 obtained using CKM unitarity 
constraints [39], we calculate /+(0) for each semileptonic mode separately and also for isospin 
averages. These are presented in Table XXIII and compared to previous measurements in 
Table XXIV [53]. The measurement of /^(O) presented here is the most precise to date. 
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TABLE XVI: Summary of results of form factor fits to the data: the normalization parameters 
/+"(0)|I4s| or /+(0)|V^d|, shape parameters (Mpoie, a, r, and ri), and the correlation coefficient p 
between the normalization and shape parameters for each fit. The last column gives the per 
degree of freedom for each fit. 



Simple pole 


/+(0)|Vc„| 


-^-^polc 


p 




per d.o.f. 






0.152(4)(1) 


1.94(4)(1) 


0.68 




1.26 




D+ nOe+Ue 


0.153(7)(5) 


1.99(10)(5) 


0.68 




0.37 




D° K~e+Ue 


0.736(7)(6) 


1 95f4l('l") 


0.78 




0.81 




D+ K°e+Ue 


0.733(11)(9) 


2.02(6)(2) 


0.78 




0.80 




D — > 7^e~^l^e 


0.152(4)(1) 


1.95(4)(2) 


0.68 




0.82 




D Ke+Ue 


0.735(7)(5) 


1.97(3)(1) 


0.78 




1.00 




Mod. pole 


f+(0)\Vcn\ 


a 


p 




per d.o.f. 






0.150(6)(1) 


0.18('12)('4) 


-0.83 




1.26 




D+ n^e+Ue 


0.151(9)(4) 


0.09(22)(12) 


-0.80 




0.35 




D° K-e+Ve 


0.733(8) (6) 


0.25f6)(2) 


-0.83 




0.94 




D+ K°e+Ue 


0.732(12)(9) 


0.12(10)(4) 


-0.82 




0.83 




D — > 7re~^i/e 


0.150(5){2) 


0.16(10)(5) 


-0.81 




0.78 




D Ke+Ue 


0.733(7)(6) 


0.21(5)(2) 


-0.83 




1.01 




ISGW2 


/+(0)|K,| 


r{ISGW2) 


P 




per d.o.f. 




D° n-e+i^e 


0.147(5)(1) 


1.98(12)(2) 


-0.80 




1.37 




D+ Tr°e+iye 


0.149(8) (4) 


1.85(22)(12) 


-0.77 




0.26 




D° -> K-e+i^e 


0.730(8)(6) 


1.56(4)(1) 


-0.81 




1.07 




D+ K°e+Ue 


0.723(12)(9) 


1.48(7)(2) 


-0.81 




0.91 




D — > Tre+i^e 


0.147(4)(2) 


1.95(10)(5) 


-0.80 




0.80 




D Ke+Ve 


0.730(7)(6) 


1.53(4)(2) 


-0.81 




1.09 




Series (2 param.) 


/+(0)|Vc,| 


ri 


1 + 1//3 - 5 




P 


per d.o.f. 


D° TT-e+I/e 


0.150(6)(1) 


-1.80(27)(5) 


1.00(11)(2) 




0.84 


1.26 


D+ n°e+Ue 


0.151(9)(4) 


-1.57(49)(26) 


0.92(17)(9) 




0.81 


0.36 


K~e+Ue 


0.734(8)(6) 


-1.96(28)(8) 


0.89(5)(2) 




0.83 


0.88 


D+ K^e+Ue 


0.733(12)(9) 


-1.40(44)(14) 


0.79(7)(3) 




0.82 


0.81 


D — > Tre+t'e 


0.151(5)(2) 


-1.75(11)(23) 


0.99(9) (4) 




0.83 


0.78 


D Ke+Ue 


0.734(7) (6) 


-1.78(24)(10) 


0.86(4) (2) 




0.83 


0.98 



TABLE XVII: Summary of the results of form factor fits for the series parametrization with three 
parameters to the data: the normalization parameters /jI^' (0)|Fcs| or /+(0)|Fcd|) the shape parame- 
ters ri, and r2, and correlation coefficients pij between the parameters determined by the fit. The 
last column gives the per degree of freedom for each fit. 



Decay 


f+mvc,\ 


ri 


r2 


1 + l/,d - 5 


poi 


P()2 


P12 


per d.o.f. 


DO -K-e+Ue 


0.152(8)(1) 


-2.0(6)(1) 


1.6(4.2)(0.8) 


0.91(24)(5) 


-0.36 


0.67 


-0.91 


1.45 


D+ TT^e+Ue 


0.144(12)(4) 


-0.3(1. 8)(1.0) 


-7.9(10.5)(5.6) 


1.36(43)(23) 


-0.46 


0.68 


-0.96 


0.29 


DO K~e+Ve 


0.745(12)(6) 


-2.4(5)(2) 


15.6(12. 8)(3.8) 


0.70(15)(5) 


-0.26 


0.71 


-0.82 


0.79 


D+ K°e+Ue 


0.744(17)(9) 


-1.9(7)(2) 


16.6(19.3)(6.2) 


0.61(21)(7) 


-0.22 


0.71 


-0.80 


0.82 


D — > Tre+fe 


0.151(7)(2) 


-1.8(6)(3) 


0.3(3.9)(1.7) 


0.98(23)(10) 


-0.38 


0.67 


-0.92 


0.84 


D Ke+i^e 


0.744(10)(6) 


-2.2(4)(2) 


16.9(11.4)(4.7) 


0.69(12)(5) 


-0.25 


0.71 


0.81 


0.93 



VII. DETERMINATION OF |Vca| AND \Vcd\ 

Using recent unquenched LQCD calculations of the form factor normalizations [26] we 
obtain \Vcq\ for each of the four semileptonic modes and for the isospin averages. These are 
presented in Table XXV for both pole models, the ISGW2 model and the series expansion 
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TABLE XVIII: Compilation of recent measurements of the ISGW2 parameters r and in units 
of GeV-i. 









Belle [14] 


2.47(15)(15) 


2.68(45) (40) 


BABAR [10] 


1.645(36) (44) 




CLEO-c (tagged) 


1.53(4)(2) 


1.95(10)(5) 



TABLE XIX: Compilation of measurements of Mpoie in D ^ Kiu(. CLEO-c (tagged) is the isospin 
averaged value; CLEO-c (untagged) is for K~e'^Ue. 



GeV/c^ 



-0.20 



Mark III [47] 
E691 [48] 
CLEO [49] 
CLEOII [50] 
E687 (Tag) [51] 
E687 (Incl) [51] 
CLEO [12] 
FOCUS [13] 
BcUc [14] 
BABAR [10] 
CLEO-c (tagged) 
CLEO-c (untagged) 



1.80 



(25) 

2.10t0-^0(20) 



2.10 



-0.20 V 
+0.40 
-0.20 



(25) 
2.00(12)(18) 



1.97; 

1.87" 



-0.43 
-0.22 
-0.11 
-0.08 



[3] 



(7) 
(7) 

1-89(5)1^^ 
1.93(5)(3) 
1.82(4)(3) 
1.884(12)(15) 
1.97(3)(1) 
1.97(3)(1) 



with two and three parameters. 

As the data do not support the physical interpretation of the shape parameter in the 
ISGW2, simple pole, and modified pole parametrizations we choose the value of |Vcg| ob- 
tained with the series expansion as our main result. Although the \ Vcq\ statistical uncertainty 
is one third larger when data is fit to the series expansion with three parameters, we choose 



TABLE XX: Compilation of measurements of Mpoie in D — Triv^. CLEO-c (tagged) is the isospin 
average value; CLEO-c (untagged) is for TT'e'^Ug. 



GeV/c^ 



CLEO (2004) [12] 
FOCUS (2004 [13]) 
Belle (2006) [14] 
CLEO-c (tagged) 



1.86 
1.91 



-1-0.10 
-0.06 
-1-0.30 
-0.15 



(5) 
(7) 



1.97(8)(4) 

1.95(4)(2) 



CLEO-c (untagged) [3] 1.87(3)(1) 
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TABLE XXI: Compilation of measurements and theoretical predictions for 
CLEO-c (tagged) are the isospin average values; CLEO-c (untagged) is for D° — 
— 7r~e+fe respectively. 



and a'^. 
K^e^Vp and 



a 



K 



a" 



FOCUS [13] 
CLEO III [12] 
Belle [14] 
BABAR [10] 
LQCD [26] 
LCSR [21] 
CQM [52] 
CLEO-c (tagged) 
CLEO-c (untagged) 



[3] 



0.28(8)(7) 
0.36(10)(5) 
0.52(8)(6) 
0.377(23)(29) 
0.50(4) 
0.07l°:J^ 

0.24 
0.21(5)(2) 
0.21(5)(3) 



0.37(25)(15) 
0.10(21)(10) 

0.44(4) 

0.30 
0.16(10)(5) 
0.37(8)(3) 



TABLE XXII: Compilation of measurements oil + 1/ fi — 5 from CLEO-c tagged and untagged [3]. 



1 + 1//3 - <5 

Ta,gg(xi Uut.agg(Hi [3] 

D^^T:~e"^Ue 0.91(24)(5) 1.30(37)(12) 
D+ TT^e+iye 1.36(43)(23) 1.58(60)(13) 
D^^K-e+Ue 0.70(15)(6) 0.62(13)(4) 
D+^K^e+i^e 0.61(21)(7) 0.51(20)(4) 
D TTe+Ue 0.98(23)(10) - 
D Ke+Ue 0.69(12)(5) 



this rather than the results obtained with the fit to two parameters to facilitate comparison 
with [3]. 

We find \Vcd\ = 0.238(12)(2)(25) for n-e+Ue, and j^dl = 0.226(21)(6)(24) 

for L)+ 7r°e+z/e. Wc find \Vcs\ = 1.020(16)(9)(106) for L)° K-e+Ue, and \Vcs\ = 
1.019(24) (13) (106) for K^e^Ve- In each case the third uncertainty in the determi- 

nation of the CKM matrix element is from theory. Averaging the and results and 
taking into account correlated and uncorrelated uncertainties we find 

\Vcd\ = 0.236 ± 0.010 ± 0.002 ± 0.025 (46) 

and 

\Vcs\ = 1.019 ± 0.013 ± 0.009 ± 0.106, (47) 

where the uncertainties are statistical, systematic and theoretical, respectively. The theoret- 
ical uncertainty dominates and is expected to be reduced soon. We compare our measure- 
ments to other determinations in Table XXVI and Table XXVII. Our determination of \Vcs\ 
is consistent with previous measurements, is in good agreement with [3], and is the most 
precise to date. Our determination of \Vcd\ is in good agreement with the result derived from 
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TABLE XXIII: Results for /+(0) obtained from fits to five form factor parametrizations. The first 
uncertainty is statistical, the second is systematic and the third uncertainty is from the relevant 
CKM matrix element. 



Mode 


Siiii]>i(* Pole 


Mod. Pole 


1SGW2 


Scries (2 param.^ 


) Series (3 paraiii.) 


D+ n^e+Ue 
DO K-e+Ue 
D+ K°e+Ue 


0.676(20) (6) (3) 0.666(25)(6)(3) 0.651(23)(6)(3) 
0.678(32)(17)(3) 0.670{39)(17){3) 0.660(35)(17)(3) 
0.756(7)(6)(0) 0.753(8)(6)(0) 0.750{8)(6)(0) 
0.753(11)(9)(0) 0.752(12)(9)(0) 0.748(12)(9)(0) 


0.667(26)(6)(3) 
0.672(40){17)(3) 
0.755(8)(6)(0) 
0.753(13)(9)(0) 


0.675(34)(6)(3) 
0.640(57)(16)(3) 
0.765(12)(7)(0) 
0.764(18)(10)(0) 


D 7if;+;/,, 
D — K< //, 


0.676(17)(7)(3) 
0.7r,6(6)(0)(0) 


0.667(21)(7)(3) 
0.7.5:i(7)(6)(0) 


0.6.'')3(19)(7)(3) 
0.7,-,0(7)((i)(0) 


0.668(21)(7)(3) 
0.7.1 i(7)((i)(0) 


0.669(29)(7)(3) 
0.7(il(10)(6)(0) 



TABLE XXIV: Compilation of measurements and theoretical predictions for f^{0) and /^(O). For 
the experimental measurements the first uncertainty is statistical, and the second is systematic. 
For this work the third uncertainty is from the relevant CKM matrix element. For BABAR the 
third uncertainty includes contributions from — K~Tr'^), tjjo and \Vcs\- 



/f(0) 



/;(o) 



LQCDl [24] 
QCD SR [21] 
LCSRl [22] 
LCSR2 [41] 
ISGW2 [20] 
LQCD2 [26] 
Belle [14] 
BABAR [10] 
CLEO-c (tagged) 
CLEO-c average 



0.66(4)(1) 

0.60(2) 
0.785(11) 
0.67(20) 

1.23 
0.73(3)(7) 
0.695(7)(22) 
0.727(7)(5)(7) 
0.764(10) (6) (0) 
0.763(7) (6) (0) 



0.57(6)(2) 
0.50(1) 
0.65(11) 
0.67(19) 

0.64(3)(6) 
0.624(20) (30) 

0.669(29) (7) (3) 
0.629(22) (7) (3) 



neutrino-nucleon scattering, it is consistent with [3] and is the most precise determination 
from D meson semileptonic decay to date. 

We also extract the ratio jVcdl/lKsl from the ratio of measured form factors. Prom the 
simultaneous quadratic z expansion fits to isospin conjugate pairs we obtain: 

j^'^jhfrl = 0-203 ± 0.009 ± 0.003, (48) 

|Vcs|/+ (0) 

where the uncertainties are statistical and systematic, respectively, and correlations have 
been taken into account. We can compare this result to calculations of /+(0) / /.^(O) to obtain 
the ratio of CKM elements. A recent light cone sum rules (LCSR) calculation obtains [56] 
fl{0)/f^{0) = 0.84 ± 0.04, from which we find 



cd\ 



0.242 ±0.011 ±0.004 ±0.012, (49) 



\Vcs. 

where the third uncertainty is from LCSR. This value is in reasonable agreement with [3] . 
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TABLE XXV: Results for IV^^I and \Vcd\ obtained from fits to five form- factor parametrizations. 
The first uncertainty is statistical, the second is systematic and the third is from the /+(0) LQCD 
prediction. 



Decay 


Simple Pole 


Mod. Pole 


ISGW2 


Series (2 param.) 


Series (3 parain.) 


D+ TfOe+I/e 

D° K-e+Ue 
D+ K°e+Ue 


0.238(7)(2)(25) 0.235(9){2)(25) 0.230(8)(2)(24) 0.235(9)(2)(25) 0.238(12)(2)(25) 
0.239(11)(6)(25) 0.236(14)(6){25) 0.233(12)(6)(24) 0.236(14)(6)(25) 0.226(20)(6)(24) 
1.008(10)(9)(105) 1.004(11)(9)(105) 1.000(11)(9)(104) 1.006(11)(8)(105) 1.020(16)(9)(106) 
1.004(15)(13)(104) 1.003(17)(13)(104) 0.997(16) (13) (104) 1.004(17)(13)(105) 1.019(24)(13)(106) 


D — » 7re+i^e 
D Ke+i^e 


0.238(6) (2) (25) 
1.007(8)(8)(105) 


0.235(7) (3) (24) 
1.004(9)(8)(105) 


0.230(7)(3)(24) 
0.999(9) (8) (104) 


0.234(8) (2) (25) 
1.006(9)(8)(105) 


0.236(10)(2)(25) 
1.019(13)(9)(106) 



TABLE XXVL Compilation of determinations of |V^s|. For the T^Kii'i) determination we use 
PDG2000 [54] , as PDG2004 [32] does not quote a value from this technique and subsequent PDG 
determinations [38, 39] include a result obtained from an earlier CLEO-c measurement, the data 
sample for which is a subset of the data used in this work. 



\Vcs\ 

1.04 ±0.16 
0.94lg:^^ ± 0.14 
1.00 ±0.05 ±0.11 
1.019 ±0.013 ±0.009 ±0.106 
1.015 ± 0.010 ± 0.011 ± 0.106 
1.018 ± 0.010 ± 0.008 ± 0.106 



r{Keiyi) PDG2000 [54] 
Charm tagged W decay [38] 
T{Kii^i)BESII [55] 
CLEO-c (tagged) 
CLEO-c (untagged) 
CLEO-c average 



VIII. CLEO-c AVERAGES 

In this section we compute average values of the measurements of branching fractions, 
form factors and \Vcs\ and \Vcd\, obtained in this work (tagged), with previous untagged 
CLEO-c measurements of the same quantities [3]. These average values represent the 
best determinations of the branching fractions, form factors, and \Vcs\ and jV^j with the 
CLEO-c 281 pb"^ data set. 

The analysis of the data, both in this work and in Ref . [3] , does not support the phys- 
ical interpretation of the shape parameter in the ISGW2, simple pole, and modified pole 



TABLE XXVIL Comparison of determinations of 



\Vcd\ 

ud cd [38] 0.22 ± 0.011 

CLEO-c (tagged) 0.236 ± 0.010 ± 0.002 ± 0.025 
CLEO-c (untagged) 0.217 ± 0.009 ± 0.004 ± 0.023 
CLEO-c average 0.222 ± 0.008 ± 0.003 ± 0.023 
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parametrization. Accordingly, both here and in Ref. [3], the values of \Vcq\ obtained with the 
series expansion with a quadratic term are chosen as the primary results. Therefore, in this 
section we present averages of \Vcq\ and the shape parameters only for the series expansion 
with a quadratic term. 

To allow external use of the set of partial branching fractions presented in this paper and 
in Ref. [3], we determine the full statistical and systematic uncertainty correlation matrices 
and present them in Appendix A. These matrices allow for simultaneous fits of the results in 
this work and in Ref. [3] to any form factor parametrization to obtain form factor parameters. 
They also allow for simultaneous fits with other experimental results. 

The two analyses use the same data set. The untagged analysis has a significantly higher 
efficiency, resulting in signal yields ~ 2.5 times greater than the tagged analysis, but also 
has larger backgrounds. Most of the signal events found by the tagged analysis are also 
found by the untagged analysis, and so the measurements produced by the two analyses are 
highly correlated. 

To compute averages we use error matrices to take into account the correlations between 
measurements made by the two techniques. The statistical covariance matrix between the 
two analyses has a 2 x 2 block form. The diagonal blocks are obtained from the untagged 
and tagged analyses, respectively. The off-diagonal blocks arise from correlations between 
the two analyses. As the covariance matrix is symmetric, only one off-diagonal block needs 
to be determined. 

The off-diagonal blocks are computed using a bootstrap [57] MC simulation, where 185 
data-sized MC samples are constructed. Each sample is created by randomly selecting events 
from the generic MC sample. Each event cannot appear more than once in a given sample. 
Each analysis nms on each bootstrap sample and the statistical correlation between the 
analyses can be measured and the off-diagonal block of the statistical covariance matrix 
computed. 

A systematic correlation matrix is constructed by taking each systematic uncertainty as 
either 100% correlated or uncorrelated between the two analyses. The complete four-block 
combined statistical and systematic correlation matrix is used to obtain the averages. A 
comparison of correlation matrices calculated by this technique to those determined by each 
analysis is made and good agreement is found. 

Consider first the determination of the average branching fractions. The untagged anal- 
ysis determines the branching fractions from the sum of the partial branching fractions in 
each bin. To treat each analysis similarly the derived branching fractions are computed 
for the tagged analysis from the partial branching fractions in each bin, using the partial 
rates in Table XII, corrected for the lifetime factors. The derived branching fractions are 
reported in Table XXVIII, and are consistent with the branching fractions found in Sec. V. 

To extract the statistical correlations between the two analyses, both analyses have run 
on the bootstrap samples. The analysis procedures in each case are identical to those used 
with data. For each semileptonic mode, and each bootstrap sample, we obtain the partial 
branching fractions in each q^ bin. We sum individual q^ bins to obtain the total branching 
fractions. These are used to calculate the statistical covariance matrices. 

Table XXIX gives the complete 8x8 statistical correlation matrix. The tagged internal 
block is diagonal as the branching fractions are uncorrelated with each other. The untagged 
internal block is obtained from the untagged analysis, where the correlations between the 
total branching fractions have been computed from those between the individual bins. 
The off-diagonal elements in the off-diagonal blocks are small, so we set them to zero. In so 
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TABLE XXVIII: The derived absolute branching fractions (in percent) obtained by summing 
the partial branching fractions in each bin, for CLEO-c tagged, untagged [3], and the average 
absolute branching fraction. 





Tagged 


Untagged 


Average 


TT e+fe 


0.308(13) (4) 


0.299(11)(8) 


0.304(11)(5) 




0.379(27)(23) 0.373(22)(13) 0.378(20)(12) 




3.60(5)(5) 


3.56(3)(9) 


3.60(3)(6) 




8.87(17)(21) 


8.53(13)(23) 


8.69(12)(19) 



doing, we neglect the statistical correlation from using common tag yields. Performing the 
averaging procedure with or without these elements produces a negligible difference in the 
final results. 

TABLE XXIX: The complete branching fraction statistical correlation matrix. Untagged inter- 
vals are in columns and tagged intervals are in rows. The modes are labeled by their final state 
hadrons. 





Untagged 
TT- K- 7r° 


Tagged 
TT- K- 7r° K° 


TT 

Untagged K~ 
^0 


1.00 -0.04 -0.02 -0.01 
1.00 0.00 -0.02 
1.00 -0.02 
1.00 


0.60 0.00 0.00 0.00 
0.00 0.49 0.00 0.00 
0.00 0.00 0.36 0.00 
0.00 0.00 0.00 0.43 


Tagged K~ 




1.00 0.00 0.00 0.00 
1.00 0.00 0.00 
1.00 0.00 
1.00 



To be conservative, all of the systematic uncertainties in the tagged analysis are taken to 
be fully correlated with the corresponding systematic uncertainties in the untagged analysis. 
The untagged analysis has additional systematic uncertainties which have no analog in the 
tagged analysis. Accordingly, these are taken to be uncorrelated between the two analyses. 
A covariance matrix for each systematic uncertainty is then constructed. A 39% correlation 
between the numbers of charged D and neutral D pairs in the untagged analysis is also 
taken into account. The complete branching fraction systematic correlation matrix is given 
in Table XXX. 

We fit the 8 branching fractions with the four-block statistical covariance matrix. The fit 
returns x^/d.o.f. = 6.7/4. We then repeat the fit with the four-block combined statistical 

and systematic correlation matrix to obtain the central values for the averages and the 
combined statistical and systematic uncertainty. This fit returns x^/d.o.f. = 2.1/4. The 
quadrature difference between the uncertainties obtained in the two fits is used to compute 
the systematic uncertainty. The results of the fits are reported in Table XXVIII. The 
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TABLE XXX: The complete branching fraction systematic correlation matrix. The modes are 
labeled by their final state hadrons. 





Untagged 






vr- K- vfO 


TT- K- VfO 


TT 


1.00 0.85 0.32 0.31 


0.48 0.46 0.08 0.23 


Untagged K~ 


1.00 0.39 0.37 


0.48 0.49 0.08 0.24 


ttO 


1.00 0.62 


0.24 0.24 0.23 0.18 


^0 


1.00 


0.38 0.38 0.10 0.59 


vr" 




1.00 0.98 0.33 0.60 


Tagged 




1.00 0.31 0.61 


tt" 




1.00 0.21 


^0 




1.00 



averaged branching fractions are more precise than those measured by either the tagged or 
untagged analysis. 

To determine the form factor parameters we perform a simultaneous fit to the (iVi + 
partial branching fractions, where A^i and N2 are the numbers of (f bins in the untagged 
and tagged analyses, respectively. To improve precision, we also simultaneously fit pairs of 
modes related by isospin. The corresponding correlation matrix has dimensions {Ni + A^2) 
X (A^i + N'i) = (56 X 56) as there are 5 bins for each mode in the untagged analysis and 
(9+7+10+10) bins in the tagged analysis. We give averaged results for the series expansion. 

For the fits to the distributions for each mode the untagged and tagged diagonal 
blocks are again taken from the analysis-determined correlation coefficients from the data 
yield fits. The bootstrap method determines the off-diagonal block statistical correlation 
matrix. As the q'^ binning differs between the two analyses we take into account both 
the correlations between intervals that overlap in the two analyses and the correlations 
between non-overlapping tagged— untagged bins. The off-diagonal block of the (56 x 56) 
statistical correlation matrix may be found in Tables A.l and A. 2. The diagonal untagged 
block is given in [3]. The systematic covariance matrix is constructed from the systematic 
uncertainties in each analysis. The off-diagonal block may be found in Tables A. 3 and A.4, 
the untagged diagonal block is given in [3], and the tagged diagonal block may be found in 
Tables A. 5 and A. 6. 

We first perform the fit with the statistical covariance matrix only, then repeat it with the 
combined statistical and systematic covariance matrix. The quadrature difference between 
the uncertainties obtained in the two fits is used to compute the systematic uncertainty. The 
central values for the averages are taken from the combined statistical and systematic fit and 
arc given in Tabic XXXI, and the parameter correlations from the fit with statistical and 
systematic uncertainties are given in Table XXXII. Finally, to improve precision, isospin 
constraints are imposed and we determine form factor parameters for D — > Ke'^u^ and 
D — > Tre'^fg from simultaneous fits to the respective isospin conjugate modes. Results and 
parameter correlations are shown in Tables XXXI and XXXIII, respectively. Using \ Vcs\ and 
\Vcd\ values constrained by CKM unitarity [39], the f^{0) and /+(0) averages are calculated 
and reported in Table XXIV. 
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TABLE XXXI: Measurements of /+(0)|Fcg|) ri and r2 in CLEO-c tagged, untagged [3] and the 
CLEO-c average. 





Tagged Untagged Average 


/+(0)|^cJ 

— > 7r~e"'"fe ri 
r2 


0.152(8)(1) 0.140(7)(3) 0.142(6)(2) 
-2.0(6)(1) -2.1(7)(3) -2.1(5)(1) 
1.6(4.2)(0.8) -1.2(4.8)(1.7) -0.6(3.5)(0.9) 


D° K'e+Ue ri 
r2 


0.745(12) (6) 0.747(9) (9) 0.745(8) (6) 

-2.4(5)(2) -2.4(4)(1) -2.4(4)(1) 
15.6(12.8)(3.8) 21(11)(2) 17.9(10.0)(1.9) 


f+mvc,\ 


0.144(12)(4) 0.138(11)(4) 0.140(9)(2) 
-0.3(1.8)(1.0) -0.2(1.5)(4) -0.2(1.3)(0.2) 
-7.9(10.5)(5.6) -9.8(9.1)(2.1) -9.4(7.5)(1.1) 




0.744(17)(9) 0.733(14)(11) 0.733(13)(9) 
-1.9(7)(2) -2.8(6)(2) -2.5(6)(1) 
16.6(19.3)(6.2) 32(18)(4) 26.8(15.4)(2.6) 


f+mvai\ 

D — ire^Vf^ ri 
r2 


0.151(7)(2) 0.142(5)(2) 
-1.8(6)(3) -1.7(5)(2) 
0.3(3.9)(1.7) -1.9(3.2)(1.1) 


/+(0)|K,l 

D — > Ke^Ve fi 
r2 


0.744(10)(6) 0.743(7)(6) 
-2.2(4)(2) -2.4(3)(1) 
16.9(11.4)(4.7) 20.8(8.3)(2.0) 



TABLE XXXIL The parameter correlation coefficients for /_|_(0)|ycg|) i^i and r2 between the tagged 
and untagged [3] CLEO-c analyses. The modes are labeled by their final state hadrons. 







n 






K- 












RO 








f+mVc<i\ ri 


r2 


/+(0)|Vc, 


n 


r2 


f+m\Vc,\ 


n 


r2 


/+(0)|Vi, 


n 


r2 




f+(0)\Vc,\ 


1.00 -0.38 


0.65 


0.11 


0.01 


0.00 


0.02 


0.00 


0.00 


0.07 


-0.01 


0.01 




ri 


1.00 


-0.93 


0.01 


0.00 


0.00 


0.01 


-0.03 


0.02 


0.01 


0.00 


0.00 




r2 




1.00 


0.00 


0.00 


0.01 


0.00 


0.01 


-0.01 


0.01 


-0.01 


0.01 










1.00 


-0.23 


0.59 


0.04 


0.00 


0.00 


0.16 


0.01 


-0.01 


K- 










1.00 


-0.83 


0.00 


-0.01 


0.01 


0.03 


-0.03 


0.03 














1.00 


0.00 


-0.01 


0.01 


-0.01 


0.03 


-0.01 
















1.00 


-0.43 


0.65 


0.03 


0.00 


0.00 


ttO 


ri 
r-i 














1.00 


-0.96 
1.00 


0.01 
-0.01 


-0.02 
0.01 


0.01 
-0.01 


Ro 


/+(0)|Vc«| 

n 


















1.00 


-0.24 
1.00 


0.62 
-0.81 
1.00 



The CLEO-c averages arc 

\Vcd\ = 0.222 ± 0.008 ± 0.003 ± 0.023 (50) 

and 

\Vcs\ = 1.018 ± 0.010 ± 0.008 ± 0.106, (51) 
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TABLE XXXIII: The parameter correlation coefficients for /_)_(0)|ycg|) I'l and r2 between the tagged 
and untagged [3] CLEO-c analysis with isospin constraints imposed. The modes are labeled by 
their final state hadrons. 





TT 






K 






f+mVc,\ n 






ri 


r2 


f+mvc,\ 


1.00 -0.42 


0.66 


0.13 


0.00 


0.01 


TT ri 


1.00 


-0.94 


0.01 


0.00 


0.00 






1.00 


0.00 


0.00 


0.01 


f+mvc,\ 






1.00 


-0.21 


0.55 


K n 








1.00 


-0.81 


r2 










1.00 



where the uncertainties arc statistical, systematic, and theoretical, respectively. We compare 
these averages to the untagged and tagged standalone determinations in Table XXVI and 
Table XXVII. In each case the average value is more precise than that obtained with either 
analysis. This is the most precise determination of \Vcs\ to date, and the most precise 
determination of \Vcd\ from D meson semileptonic decay to date. 

IX. SUMMARY 

In this paper we have presented precise measurements of the absolute branching fractions 
of D'^ decays to K'e^Ug and ir^e^Ug and decays to K^e^u^ and Tc'^e^u^, that agree well 
with world averages [32] and Ref . [3] . We have combined these measurements to demonstrate 
that D meson exclusive semileptonic decays to pseudoscalar final states are consistent with 
isospin invariance. 

From the spectrum of all four decay modes studied, we have made the most precise 
determinations of and r^. Our measurement of is over Aa smaller than Belle [14], and 
1.6(7 smaller than BABAR [10]. Our measurement of M^jg (a^) is in excellent agreement 
with Ref. [3], but is 2.3a larger (2. 60" smaller) than the BABAR [10] measurement, and 
2.5(7 larger (2.7o" smaller) than the Belle [14] measurement. Our determinations of r'^, M^^^^ 
and are in reasonable agreement with Ref [3] and other previous measurements. Our 
measurement of {a'^) is more than Aa (2.4(j) smaller than the LQCD fit. However, the 
discrepancy with LCQD is difficult to quantify because the covariance matrix for the LQCD 
form factors is lost during the chiral extrapolation procedure for the published analysis [26] . 
The results reported here for the series expansion parameters ri and r2 are in excellent 
agreement with Ref. [3] and agree with BABAR [10] to better than 2a with the precise level 
depending on the correlation coefficient for the BABAR ri and r2 parameters. 

Our data, and other recent measurements (e.g., Refs. [10, 12-14]) do not support the 
physical interpretation of the shape parameter in the ISGW2, simple pole, and modified 
pole parametrizations. Accordingly, the /+(0)|Vcg| values obtained when the data is fit 
with the quadratic series expansion were selected as our primary normalization results. We 
combined these values with the unitarity of the CKM matrix to make a precise determination 
of the form factor absolute magnitude f^{0) and the most precise determination of /+(0). 
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Using unqucnched LQCD predictions for f^{0) and /+(0) we have made the most precise 
determination of \Vcs\, and the most precise determination of \Vcd\ from D semileptonic 
decays to date. The results agree well with previous measurements using semileptonic decays 
including Ref [3], and agree well with charm-tagged W decay measurements of \Vcs\ and 
neutrino based determinations of \Vcd\- 

To allow external use of the set of partial branching fractions presented in this paper and 
in Ref. [3] we determined the full statistical and systematic uncertainty correlation matrices. 
These matrices allow for simultaneous fits of the results of this work and Ref [3] to any form 
factor parametrization to obtain form factor parameters. They also allow for simultaneous 
fits with other experimental results. 

Finally, we averaged values of the measurements obtained in this work with Ref. [3]. 
These averages represent the best determinations of the branching fractions, form factors 
and \Vcs\ and \Vcd\ with the CLEO-c 281 pb~^ data set [58]. They are the most precise 
measurements of the absolute branching fractions of decays to K^e^u^, and Tr^e+z/g and 
decays to K'^e^v^ and 7r°e+z/e, and the most precise direct determination of \Vcs\ and 
the most precise determination of \Vcd\ from D semileptonic decay. 

CESR has recently collected a larger t/'(3770) data sample. It is expected that this sample 
will result in a further improvement in measurements of and semileptonic branching 
fractions, measurements of the decay form factors, and the CKM matrix elements \Vcs\ and 
\Vcd\ [59]. 
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APPENDIX A: CORRELATION MATRICES 

To allow external use of the set of partial branching fractions presented in this paper 
and [3] this Appendix contains the statistical and systematic uncertainty correlation ma- 
trices. These matrices allow for simultaneous fits of these results with other experimental 
results to obtain form factor parameters. Sec. VIII describes the procedures that have been 
used to obtain these matrices. 



43 



TABLE A.l: The untagged-tagged block of the statistical correlation matrix obtained from the 
bootstrap procedure . Untagged intervals are in columns and tagged intervals are in rows. 
The lines indicate mode boundaries. The modes are labeled by their final state hadrons. Within 
each submode, the intervals are ordered from lowest to highest (part I). 





TT 


Tagged 

K- 




0.46 0.20 -0.01 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.35 0.30 -0.02 0.00 0.00 0.00 0.00 0.00 
0.00 -0.02 0.28 0.53 0.01 -0.01 0.00 0.00 0.00 
0.00 0.00 -0.01 -0.02 0.43 0.20 0.00 0.00 -0.01 
0.00 0.00 0.00 0.00 -0.01 0.35 0.26 0.15 0.35 


-0.01 -0.02 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 


K- 

Untagged 


-0.03 -0.01 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 


0.25 0.27 0.04 -0.02 -0.01 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.44 0.37 0.08 -0.03 -0.01 0.00 0.00 0.00 
0.00 0.00 -0.03 -0.03 0.28 0.48 0.13 -0.03 -0.01 0.00 

0.00 0.00 0.00 0.00 -0.02 -0.03 0.28 0.47 0.17 -0.02 
0.00 0.00 0.00 0.00 0.00 0.00 -0.03 -0.04 0.25 0.23 




0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 -0.04 -0.03 -0.02 -0.04 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 




0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 -0.02 -0.01 -0.01 -0.02 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -0.02 -0.01 



TABLE A. 2: The untagged-tagged block of the statistical correlation matrix obtained from the 
bootstrap procedure . Untagged intervals are in columns and tagged intervals are in rows. 
The lines indicate mode boundaries. The modes are labeled by their final state hadrons. Within 
each submode, the intervals are ordered from lowest to highest (part II). 







Tagged 


TT 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 -0.03 -0.04 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -0.02 -0.01 


K- 

Untagged 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -0.02 -0.02 




0.28 0.18 -0.02 0.00 0.00 0.00 0.00 
0.00 0.31 0.27 -0.03 0.00 0.00 0.00 
0.00 -0.03 0.12 0.36 0.00 -0.01 0.00 
0.00 0.00 -0.01 -0.03 0.41 0.13 -0.02 

0.00 0.00 0.00 0.00 -0.03 0.22 0.30 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 

0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -0.02 -0.01 




0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 0.00 0.00 -0.01 -0.02 


0.34 0.28 0.04 -0.01 -0.01 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.26 0.21 0.09 -0.02 -0.01 0.00 0.00 0.00 
0.00 0.00 -0.02 -0.01 0.22 0.39 0.20 -0.02 -0.01 0.00 
0.00 0.00 0.00 0.00 -0.02 -0.03 0.23 0.26 0.17 -0.03 
0.00 0.00 0.00 0.00 0.00 0.00 -0.02 -0.03 0.33 0.27 
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TABLE A. 3: The untagged-tagged block of the systematic correlation matrix. Untagged inter- 
vals are in columns and tagged intervals are in rows. The lines indicate mode boundaries. The 
modes are labeled by their final state hadrons. Within each submode, the intervals are ordered 
from lowest to highest (part I). 





1 

7r~ 


^agged 

K- 




0.34 0.34 0.33 0.33 0.33 0.32 0.32 0.31 0.29 

0.44 0.47 0.47 0.47 0.47 0.47 0.47 0.44 0.44 
0.39 0.40 0.39 0.39 0.38 0.38 0.38 0.36 0.35 
0.31 0.32 0.31 0.31 0.31 0.31 0.31 0.29 0.28 
0.26 0.26 0.26 0.26 0.25 0.25 0.25 0.24 0.23 


0.29 0.32 0.33 0.33 0.33 0.33 0.32 0.31 0.30 0.21 

0.32 0.39 0.42 0.44 0.45 0.45 0.45 0.44 0.44 0.28 
0.33 0.36 0.37 0.38 0.38 0.38 0.38 0.36 0.36 0.24 
0.25 0.28 0.30 0.31 0.31 0.31 0.31 0.30 0.29 0.20 
0.21 0.24 0.25 0.25 0.25 0.25 0.25 0.24 0.24 0.16 


Untagged 


0.37 0.38 0.37 0.37 0.37 0..36 0.36 0.35 0.33 
0.33 0.34 0.33 0.33 0.33 0.32 0.32 0.31 0.29 
0.33 0.33 0.33 0.33 0.32 0.32 0.32 0.31 0.29 
0.27 0.27 0.26 0.26 0.26 0.25 0.25 0.25 0.23 
0.26 0.27 0.27 0.27 0.26 0.26 0.26 0.25 0.23 


0.37 0.39 0.39 0.39 0.40 0.39 0.39 0..38 0.37 0.27 
0.33 0.35 0.35 0.35 0.35 0.35 0.34 0.34 0.33 0.24 
0.34 0.35 0.35 0.35 0.35 0.35 0.34 0.34 0.32 0.24 
0.28 0.29 0.29 0.29 0.29 0.28 0.28 0.27 0.26 0.20 
0.26 0.27 0.28 0.28 0.28 0.28 0.27 0.27 0.26 0.19 




0.19 0.18 0.18 0.17 0.17 0.16 0.16 0.16 0.14 
0.14 0.15 0.16 0.16 0.16 0.15 0.15 0.15 0.14 
0.09 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.09 
0.23 0.23 0.23 0.23 0.22 0.22 0.22 0.21 0.20 
0.10 0.11 0.11 0.12 0.11 0.12 0.12 0.11 0.11 


0.21 0.21 0.20 0.20 0.19 0.19 0.18 0.17 0.16 0.12 
0.14 0.15 0.16 0.16 0.16 0.16 0.16 0.14 0.14 0.09 
0.08 0.09 0.10 0.10 0.10 0.10 0.10 0.09 0.09 0.06 
0.19 0.21 0.21 0.21 0.22 0.22 0.21 0.20 0.20 0.14 
0.07 0.09 0.10 0.11 0.11 0.11 0.11 0.11 0.12 0.08 




0.24 0.24 0.23 0.23 0.22 0.21 0.21 0.21 0.19 
0.29 0.29 0.28 0.28 0.28 0.27 0.27 0.26 0.24 
0.30 0.30 0.30 0.30 0.29 0.29 0.29 0.28 0.26 
0.31 0.31 0.31 0.31 0.30 0.30 0.30 0.29 0.27 
0.23 0.24 0.23 0.23 0.23 0.22 0.22 0.21 0.20 


0.27 0.27 0.26 0.26 0.26 0.25 0.25 0.23 0.22 0.17 
0.29 0.31 0.30 0.30 0.30 0.30 0.29 0.28 0.27 0.20 
0.29 0.31 0.32 0.32 0.32 0.31 0.31 0.30 0.29 0.21 
0.30 0.32 0.32 0.33 0.33 0.32 0.32 0.31 0.30 0.21 
0.22 0.24 0.24 0.24 0.24 0.24 0.23 0.22 0.22 0.15 



TABLE A. 4: The untagged-tagged block of the systematic correlation matrix. Untagged inter- 
vals are in columns and tagged intervals are in rows. The lines indicate mode boundaries. The 
modes are labeled by their final state hadrons. Within each submode, the q^ intervals are ordered 
from lowest to highest (part II). 







Tagged 




0.05 0.07 0.08 0.10 0.11 0.13 0.16 
0.05 0.09 0.12 0.14 0.16 0.20 0.25 
0.05 0.08 0.10 0.12 0.13 0.16 0.19 
0.03 0.05 0.07 0.08 0.10 0.12 0.15 
0.02 0.04 0.05 0.07 0.08 0.10 0.12 


0.16 0.19 0.20 0.20 0.21 0.21 0.21 0.22 0.20 0.17 
0.14 0.22 0.24 0.26 0.27 0.28 0.28 0.31 0.28 0.23 
0.16 0.21 0.22 0.22 0.23 0.24 0.24 0.25 0.22 0.19 
0.11 0.15 0.17 0.17 0.18 0.19 0.19 0.20 0.18 0.15 
0.09 0.12 0.13 0.13 0.14 0.14 0.14 0.15 0.14 0.11 


K- 

Untagged 


0.05 0.07 0.09 0.11 0.12 0.15 0.18 
0.05 0.07 0.08 0.10 0.11 0.13 0.16 
0.05 0.07 0.08 0.09 0.11 0.13 0.16 
0.04 0.05 0.06 0.07 0.08 0.10 0.12 
0.03 0.05 0.06 0.07 0.08 0.10 0.13 


0.17 0.21 0.21 0.22 0.23 0.23 0.23 0.23 0.21 0.18 
0.15 0.18 0.19 0.20 0.20 0.21 0.21 0.21 0.19 0.16 
0.16 0.19 0.20 0.20 0.21 0.21 0.21 0.21 0.19 0.16 
0.13 0.15 0.16 0.16 0.16 0.17 0.17 0.17 0.15 0.13 
0.09 0.12 0.12 0.13 0.13 0.14 0.14 0.14 0.13 0.11 




0.23 0.23 0.24 0.24 0.23 0.22 0.20 
0.16 0.17 0.18 0.18 0.18 0.18 0.18 
0.17 0.18 0.18 0.18 0.18 0.17 0.15 
0.32 0.33 0.33 0.34 0.33 0.32 0.29 
0.24 0.25 0.25 0.25 0.24 0.23 0.19 


0.18 0.18 0.17 0.17 0.17 0.17 0.17 0.16 0.15 0.13 
0.13 0.15 0.15 0.15 0.15 0.16 0.15 0.15 0.14 0.12 
0.07 0.09 0.09 0.09 0.09 0.10 0.09 0.10 0.09 0.08 
0.13 0.15 0.15 0.16 0.16 0.17 0.17 0.17 0.15 0.13 
0.04 0.06 0.07 0.08 0.08 0.08 0.08 0.09 0.08 0.07 




0.07 0.08 0.09 0.11 0.12 0.14 0.16 

0.07 0.09 0.11 0.12 0.14 0.16 0.20 
0.07 0.09 0.11 0.13 0.15 0.17 0.20 
0.08 0.10 0.12 0.13 0.15 0.18 0.21 
0.06 0.07 0.09 0.10 0.12 0.14 0.17 


0.52 0.52 0.51 0.51 0.50 0.50 0.49 0.47 0.46 0.39 

0.51 0.53 0.52 0.52 0.52 0.52 0.51 0.49 0.48 0.40 
0.50 0.51 0.51 0.52 0.51 0.51 0.51 0.49 0.47 0.40 
0.46 0.49 0.49 0.49 0.49 0.49 0.48 0.47 0.45 0.38 
0.34 0.36 0.36 0.36 0.36 0.36 0.35 0.34 0.33 0.28 
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TABLE A. 5: The tagged block of the systematic correlation matrix. The lines indicate mode 
boundaries. The modes are labeled by their final state hadrons. Within each submode, the 
intervals are ordered from lowest to highest (part I). 



Tagged 



Tagged 



1.00 0.93 0.91 0.90 0.89 0.88 0.88 0.85 0.80 
1.00 0.94 0.94 0.92 0.92 0.92 0.88 0.84 
1.00 0.93 0.92 0.91 0.91 0.87 0.83 
1.00 0.92 0.92 0.92 0.88 0.84 
1.00 0.90 0.90 0.86 0.83 
1.00 0.90 0.86 0.83 
1.00 0.86 0.83 
1.00 0.79 
1.00 



0.85 
0.81 
0.79 
0.77 
0.76 
0.75 
0.75 
0.73 
0.66 



0.91 0, 
0.89 0, 
0.88 0, 
0.87 0, 
0.85 0, 
0.84 0, 
0.84 0. 
0.82 0. 
0.75 0. 



92 0.92 
92 0.93 
91 0.92 
91 0.92 
89 0.90 

88 0.90 

89 0.90 
86 0.87 
80 0.81 



K' 
0.92 0, 
0.94 0, 
0.93 0, 
0.93 0, 
0.91 0, 
0.91 0, 
0.91 0. 
0.88 0. 
0.82 0. 



91 0.90 
93 0.92 
93 0.91 
93 0.91 
91 0.90 
91 0.90 
91 0.90 
87 0.86 
83 0.82 



0.87 0, 
0.90 0, 
0.89 0, 
0.89 0, 
0.88 0, 
0.88 0, 
0.88 0. 
0.84 0. 
0.80 0. 



84 0.60 
87 0.61 
87 0.60 
87 0.60 
86 0.59 
86 0.58 
86 0.58 
82 0.56 
79 0.53 



1.00 



0.94 0. 
1.00 0. 
1. 



91 0.90 
96 0.95 
00 0.97 
1.00 



K- 



0.88 0. 
0.94 0. 
0.96 0. 
0.97 0. 
1.00 0. 
1. 



86 0.85 
92 0.91 

95 0.94 

96 0.94 
96 0.95 
00 0.94 

1.00 



0.82 
0.89 
0.92 
0.92 
0.93 
0.92 
0.91 
1.00 
1 



.77 0.62 
,84 0.64 
,88 0.65 
,89 0.65 
,90 0.65 
,89 0.64 
,88 0.63 
,87 0.63 
,00 0.60 
1.00 
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TABLE A. 6: The tagged block of the systematic correlation matrix. The lines indicate mode 
boundaries. The modes are labeled by their final state hadrons. Within each submode, the 
intervals are ordered from lowest to highest (part II). 



0.19 
0.18 
0.18 
0.17 
0.17 
0.17 
0.17 
0.16 
0.15 



0.25 0.30 
0.25 0.30 
0.24 0.30 
0.24 0.30 
0.24 0.29 
0.23 0.29 
0.23 0.29 
0.23 0.28 
0.21 0.26 



JT^ 

0.34 0.40 
0.35 0.40 
0.34 0.40 
0.34 0.40 
0.34 0.40 
0.34 0.39 
0.34 0.39 
0.32 0.38 
0.31 0.36 



Tagged 



0.46 
0.48 
0.47 
0.47 
0.47 
0.47 
0.47 
0.45 
0.43 



0.56 
0.58 
0.58 
0.58 
0.57 
0.57 
0.57 
0.55 
0.52 



0.46 
0.43 
0.41 
0.40 
0.40 
0.38 
0.38 
0.38 
0.33 



0.54 0.55 
0.53 0.55 
0.51 0.54 
0.51 0.54 
0.50 0.53 
0.49 0.53 
0.49 0.53 
0.48 0.51 
0.44 0.47 



0.57 0. 
0.57 0. 
0.56 0. 
0.56 0. 
0.55 0. 
0.55 0. 
0.55 0. 
0.53 0. 
0.50 0. 



58 0.59 

59 0.60 
58 0.60 
58 0.60 
57 0.59 
56 0.58 
56 0.58 
54 0.56 
51 0.53 



0.58 0.59 
0.60 0.62 
0.59 0.62 
0.59 0.62 
0.58 0.61 
0.58 0.61 
0.58 0.61 
0.56 0.58 
0.53 0.56 



0.55 0.46 

0.57 0.48 

0.56 0.47 

0.56 0.47 

0.56 0.47 

0.55 0.47 

0.55 0.47 

0.53 0.45 

0.51 0.43 








19 





23 





27 


0.31 


0.35 





40 


0.48 


0.53 





56 





55 





55 





55 


0.55 


0.55 





52 





50 


0.42 







19 





24 





29 


0.33 


0.38 





44 


0.53 


0.52 





57 





58 





59 





59 


0.60 


0.59 





58 





55 


0.46 







19 





25 





30 


0.34 


0.39 





46 


0.56 


0.49 





56 





58 





59 





60 


0.61 


0.61 





61 





57 


0.48 







19 





25 





30 


0.34 


0.40 





46 


0.57 


0.48 





56 





58 





59 





60 


0.62 


0.61 





62 





58 


0.48 


K- 





18 





24 





30 


0.34 


0.40 





47 


0.57 


0.47 





55 





58 





59 





60 


0.62 


0.61 





63 





58 


0.49 







18 





24 





29 


0.34 


0.40 





47 


0.57 


0.45 





54 





57 





59 





60 


0.61 


0.61 





63 





58 


0.48 







18 





24 





29 


0.34 


0.39 





46 


0.56 


0.44 





54 





56 





58 





59 


0.61 


0.60 





62 





57 


0.48 







17 





23 





28 


0.33 


0.38 





45 


0.55 


0.42 





52 





54 





56 





57 


0.59 


0.58 





60 





55 


0.47 







16 





22 





27 


0.32 


0.37 





44 


0.54 


0.39 





49 





52 





54 





55 


0.57 


0.57 





59 





54 


0.46 







12 





16 





19 


0.22 


0.26 





30 


0.37 


0.32 





37 





38 





39 





39 


0.40 


0.40 





40 





38 


0.32 



Tagged 



1.00 



0.98 0.97 
1.00 0.98 
1.00 



0.95 0.91 
0.96 0.93 
0.97 0.95 
1.00 0.95 
1.00 



0.84 
0.87 
0.90 
0.91 
0.93 
1.00 



0.67 
0.72 
0.77 
0.80 
0.85 
0.89 
1.00 



0.15 
0.17 
0.19 
0.22 
0.24 
0.28 
0.33 



0.15 0.15 
0.19 0.19 
0.22 0.22 
0.25 0.26 
0.28 0.29 
0.33 0.34 
0.39 0.41 



0.15 0.15 0.15 

0.19 0.19 0.19 

0.23 0.23 0.23 

0.26 0.26 0.27 

0.30 0.30 0.31 

0.35 0.35 0.36 

0.42 0.43 0.44 



0.15 0.14 
0.19 0.19 
0.23 0.23 
0.27 0.27 
0.31 0.31 
0.36 0.37 
0.44 0.45 



0.14 0.12 
0.18 0.15 
0.22 0.18 
0.25 0.21 
0.29 0.25 
0.34 0.29 
0.42 0.35 



1.00 



0.96 0.94 
1.00 0.96 
1.00 



0.93 0.91 0.90 
0.96 0.95 0.94 
0.96 0.95 0.95 
1.00 0.95 0.95 
1.00 0.95 
1.00 



0.89 0.84 
0.93 0.89 
0.94 0.91 
0.94 0.92 
0.94 0.91 
0.94 0.92 
1.00 0.91 
1.00 



0.82 0.69 
0.87 0.73 
0.88 0.74 
0.88 0.74 
0.88 0.74 
0.88 0.74 
0.88 0.74 
0.86 0.72 
1.00 0.69 
1.00 
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